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ABSTRACT 

Space  Shuttle  Orbiter  landings  indicate  both  long  term  directional  instability  and  the 
potential  for  pilot  induced  oscillations  during  landing  and  rollout  before  nosewheel  touchdown. 
The  Orbiter' s  Flight  Control  System  requires  improvements  to  increase  directional  control  in  the 
two  point  stance  (after  main  gear  touchdown  with  the  nose  in  the  air).  A  number  of  modifications 
are  proposed  to  improve  directional  control.  This  thesis  describes  the  control  deficiency,  potential 
improvements  to  the  Flight  Control  System  (FCS),  and  evaluates  a  number  of  these 
improvements.  The  evaluation  was  performed  by  modeling  the  Orbiter' s  postlanding 
lateral/directional  control  laws  using  a  commercially  available  engineering  software  package 
known  as  MATLAB  5.0.  Directional  control  of  the  Orbiter  was  evaluated  with  and  without  the 
proposed  modification  to  obtain  a  comparison  of  control  response. 

Initial  evaluation  of  future  Orbiter  FCS  modifications  could  be  performed  using 
commercially  available  engineering  software  packages  such  as  MATLAB;  rather  than  costly  full- 
up  Orbiter  simulators.     A  low  cost  initial  evaluation  of  changes  may  save  NASA  resources. 
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I.    INTRODUCTION 

The  Space  Shuttle,  operated  by  the  National  Aeronautics  and  Space  Administration 
(NASA),  is  the  only  manned  space  flight  vehicle  currently  operated  by  the  United  States.  The 
Space  Shuttle  consists  of  the  Orbiter,  External  Tank  (ET),  and  two  Solid  Rocket  Boosters 
(SRB's).  The  Orbiter  contains  the  astronauts  and  payloads.  It  is  the  only  piece  of  the  Space 
Shuttle  that  completes  the  entire  mission.  The  Orbiter  returns  from  space  to  land  similar  to  a 
conventional  gliding  airplane.  The  continued  safe  operation  of  the  Space  Shuttle  is  the  number 
one  concern  of  the  Space  Shuttle  Program  (SSP).  A  recent  safety  concern  is  the  directional 
control  of  the  Orbiter  during  one  segment  of  the  landing/rollout  phase  [Ref  1].  This  perceived 
directional  control  deficiency  is  described  below. 

A.    BACKGROUND 

During  the  landing  and  rollout  phase  of  flight,  the  gains  associated  with  the  piloting  task 
are  greatly  increased.  The  margin  for  error  is  vastly  decreased  as  pitch,  roll,  and  yaw  requirements 
become  more  stringent  through  the  flare-to-touchdown  and  subsequent  rollout.  During  this 
critical  phase,  pilot  control  commands  will  be  more  frequent  and  often  of  greater  magnitude  to 
make  timely  corrections  to  flight  path  and  counter  external  perturbations  such  as  crosswinds,  wind 
shear,  or  gusty  conditions.  These  increased  demands  on  the  control  system  can  occasionally  result 
in  a  condition  known  as  Pilot  Induced  Oscillation  (PIO). 

PIO  occurs  when  a  pilot  attempts  a  particular  control  task  which  is  beyond  the  capabilities 
of  the  aircraft  control  system.  PIO  can  only  occur  with  the  pilot  actively  in  the  vehicle  control 
loop  and  is  therefore  not  a  stability  concern.  A  PIO  is  typically  manifested  by  an  uncommanded 
oscillatory  motion  that  can  only  be  stopped  when  the  pilot  reduces  control  inputs  or  'backs  out  of 
the  loop.'  The  severity  of  PIO  can  range  from  an  annoyance  to  complete  loss  of  control.  The 
result  of  a  PIO  situation  can  often  be  unpredictable.  What  might  be  an  annoyance  during  landing 
one  day,  might  be  disastrous  on  the  next.  PIO  can  result  from  many  different  root  causes 
including  an  excess  amount  of  friction  in  a  mechanical  control  system,  actuator  rate  limiting,  or 
transport  delays  in  today's  most  advanced  flight  control  systems.  In  all  cases,  a  large  phase  lag 
from  the  control  input  to  aircraft  response,  lies  at  the  heart  of  the  problem. 

The  Orbiter  incorporates  a  'fly-by-wire'  control  system  to  command  the  flight  control 
surfaces,  depicted  in  Figure  1 . 


Inboard  Elevon 


Body  Flap 


Rudder/Speedbrake 


Outboard  Elevon 


Figure  1 .  Orbiter  Flight  Control  Surfaces 


The  flight  controls,  control  stick,  rudder  pedals,  and  speed  brake  switch  are  not 
mechanically  connected  to  their  respective  control  surfaces.  Orbiter  control  is  afforded  through  a 
myriad  of  feedback  control  loops,  filters,  gain  multipliers,  integrators,  differentiators,  etc.  The 
goal  of  these  components,  in  conjunction  with  a  feel  system,  is  to  supply  the  pilot  with  an  Orbiter 
which  responds  acceptably  to  all  control  inputs.  Fly-by-wire  control  systems,  if  not  carefully 
designed,  can  be  susceptible  to  PIO. 

The  directional  control  problem  of  interest  is  most  evident  after  touchdown  while  the 
Orbiter  is  in  the  two  point  stance  with  the  mainwheels  on  the  ground  and  the  nosewheel  in  the  air. 
This  flight  segment,  from  initial  touchdown  until  denotation  is  complete,  is  approximately  20 
seconds  in  duration.  With  dispersed  touchdown  speeds  between  185  and  225  knots,  a  small 
deviation  in  track  can  result  in  a  runway  departure  very  quickly.  Recent  Orbiter  landings  have 
shown  that  there  is  a  directional  PIO  problem  in  the  two  point  stance[Ref  1].  The  problem  seems 
to  abate  after  nose  gear  touchdown  due  to  the  additional  directional  control  afforded  by  the  Nose 
Wheel  Steering  (NWS)  system.  The  PIO  is  associated  with  the  task  of  capturing  and  tracking  the 
runway  centerline  and  may  be  exaggerated  in  conditions  of  greater  crosswind.  The  possibility  of 
the  Orbiter  becoming  directionally  uncontrollable  and  departing  the  runway  is  unacceptable.  This 
problem  raises  some  critical  questions  including: 


1 .  Are  there  operational  impacts?  Must  the  crosswind  limit  be  reduced  to  protect  against 
PIO?  Should  Shuttle  landings  require  wider  runways  and/or  lakebed  landings? 

2.  Are  there  hardware  modifications  such  as  changes  to  the  rudder  actuators  that  would 
improve  directional  control? 

3.  Are  there  potential  software  modifications  to  the  control  laws  which  will  reduce  or 
eliminate  this  PIO  condition? 

4.  If  PIO  is  evident  in  a  nominal  landing,  what  would  happen  in  the  case  of  a  blown 
main  gear  tire? 

1.  Operational  Impacts 

Reducing  the  crosswind  limit  may  reduce  the  potential  for  direction  PIO;  however,  it  will 
not  eliminate  the  problem.  In  fact,  the  landing  which  experienced  the  most  significant  directional 
PIO  had  a  2  knot  crosswind  [Ref  2].  Other  flight  data  also  shows  that  even  with  minimal 
crosswind,  directional  PIO  can  occur. 

In  addition,  there  is  appreciable  pressure  to  increase  the  crosswind  limits  for  the  following 
reason.  The  commitment  of  the  United  States  to  building  the  International  Space  Station  has 
extended  the  mission  of  the  SSP  from  pure  science  to  upholding  national  political  obligation.  The 
proposed  International  Space  Station  will  be  in  a  high  inclination  orbit  which  will  reduce  launch 
opportunities  and  decrease  the  length  of  launch  windows.  Landing  weather  conditions  are  a 
concern  on  launch  day  since  the  potential  for  an  emergency  landing  such  as  a  Return  To  Launch 
Site  (RTLS)  or  Trans-Atlantic  Landing  (TAL)  can  occur  on  launch  day.  Meeting  the  Space 
Station  construction  schedule  would  be  much  easier  if  launches  and  landings  could  be  performed 
in  a  greater  variety  of  wind  conditions. 

Requiring  lakebed  landings  for  all  Shuttle  flights  is  unacceptable.  Orbiter  landing  mass 
properties  are  greatly  restricted  for  lakebed  landings  to  maintain  nose  gear  loads  within  limits 
[Ref.  3].  Most  Orbiter  landing  mass  properties  require  concrete  runways.  Further,  landing  on  a 
lakebed  requires  landing  at  either  Edwards  Air  Force  Base  (EAFB)  in  California  or  Northrup, 
New  Mexico  which  is  impossible  for  most  ascent  aborts  modes. 

The  Shuttle  Landing  Facility  (SLF)  at  Kennedy  Space  Center  (KSC)  and  the  concrete 
runway  used  for  Shuttle  landings  at  EAFB  are  300  ft.  wide.  The  SLF  has  50  ft.  of  usable  shoulder 
on  each  side.  In  comparison,  some  contingency  and  abort  runways  are  less  than  150  ft.  wide. 
Obviously,  it  is  impossible  for  the  SSP  to  widen  all  runways  available  for  landing.   Restricting  the 


Orbiter  from  landing  on  runways  considered  too  narrow,  greatly  reduces  abort  coverage  (no  where 
to  landing  with  an  SSME  failed)  to  an  unacceptable  level. 

2.  Hardware  Modifications 

Hardware  modifications,  such  as  redesign  of  the  rudder  actuator  to  increase  deflection 
rate,  may  reduce  susceptibility  to  PIO.  However,  hardware  modifications  are  expensive,  require 
recertification,  and  demand  vehicle  down  time  to  implement  the  changes. 

3.  Software  Modifications 

The  obvious  long  term  solution  to  the  PIO  problem  lies  in  redesigning  the  existing  flight 
control  laws  to  afford  better  directional  control  in  a  wider  variety  of  wind  conditions. 

A  redesign  of  this  magnitude  would  traditionally  require  the  control  engineer  to  estimate 
the  new  control  laws  followed  by  experimental  Verification  and  Validation  (V&V)  of  the  new 
control  laws  in  the  Vertical  Motion  Simulator  (VMS)  located  at  the  NASA  Ames  Research  Center 
(ARC)  located  in  Sunnyvale,  California.  This  procedure  typically  entails  qualified  astronauts 
traveling  from  the  Johnson  Space  Center  (JSC)  in  Houston,  Texas  to  the  ARC  VMS  to  fly 
numerous  simulated  approaches  in  varying  wind  conditions.  This  evolution  can  be  time 
consuming  and  costly. 

With  the  advent  of  faster  more  powerful  personal  computers  and  workstations,  and  the 
availability  of  very  capable  software  packages  which  can  operate  in  both  the  frequency  and  time 
domains,  the  possibility  of  writing  and  evaluating  new  control  laws  immediately  has  become  a 
reality.  This  breakthrough  can  eliminate  the  costly  iterative  process  of  tweaking  and  reevaluating 
control  laws  in  an  expensive  motion  base  simulator.  The  final  V&V  is  still  necessary,  but  the  total 
time  and  cost  are  drastically  reduced. 

The  use  of  simulation  to  evaluate  control  laws  rely  heavily  on  the  numerous  models  used 
to  develop  the  simulation.  A  normal  flight  test  program  would  not  rely  solely  on  simulation 
results  for  solutions  to  PIO  problems.  The  Shuttle  Program  must  rely  on  these  simulations  results 
due  to  the  nature,  cost,  and  risk  of  flying  Shuttle  flights  only  to  evaluate  control  law  changes. 
Flight  data  is  gathered  from  each  flight;  however,  flight  conditions  do  not  test  the  extreme  ends  of 
the  envelope.  There  is  no  possibility  of  performing  a  complete  flight  test  program. 


B.         PURPOSE 


This  thesis  has  two  objectives: 

(1)  Develop  changes  to  the  control  laws  to  reduce  the  directional  PIO  tendency  present 
between  mainwheel  touchdown  and  completion  of  derotation  while  the  Orbiter  is  in 
the  two  point  attitude. 

(2)  Evaluate  these  changes  using  a  Commercial  Off-The-Shelf  (COTS)  PC-based 
software  package. 


II.  ORBITER  DIRECTIONAL  STABILITY  AND  CONTROL 

The  nominal  Obiter  rollout  profile,  directional  stability,  and  directional  control  in  the  two 
point  stance  are  described  in  this  chapter. 

A.         NOMINAL  ROLLOUT  PROFILE 

The  nominal  touchdown  and  rollout  profile,  with  approximate  speeds  and  time  from 
touchdown,  is  described  in  Table  1  [Ref.  2,  3  and  4]. 

Table  1 .  Nominal  Touchdown/Rollout  Profile 


Event 

Airspeed  (KEAS) 

Time  (sec.) 

Touchdown* !) 

205(195) 

0 

Brakes  functional*  ' 

200 

1.9 

Drag  chute  deploy(3) 

195 

3.0 

Commence  auto  denotation'4' 

185 

11.0 

Commence  manual 
derotation(5) 

175 

16.0 

Drag  chute  fully  disreefed 

175 

16.0 

Nosewheel  touchdown 

145  to  155 

19.0 

First  brake  application'6' 

140  (KGS) 

23.0 

Drag  chute  Jettison(7) 

60  ±  20  (KGS) 

34.0 

Wheels  stop 

0 

50.0 

Notes:    1 .    Targeted  touchdown  speed  of  205  KEAS  for  Orbiter  gross  weights  > 
220K  lbs.  (195  KEAS  for  Orbiter  weights  <  220K  lbs). 

2.  Built  In  Test  (BIT)  requires   1.9  sec.  after  Weight  On  Wheels  Lock-On 
(WOWLON). 

3.  Chute  deploy  after  main  gear  touchdown. 

4.  Nominal  derotation  is  automatic  with  beep  trim  activation  targeting   1.5 
deg/sec. 

5.  Only  performed  if  beep  trim  fails. 

6.  Brake  initiation  at  140  KGS  or  5,000  ft.  runway  remaining,  whichever  occurs 
first. 


7.    If  drag  chute  is  not  jettisoned  by  40  KGS,  the  chute  will  remain  on  to 
minimize  damage  to  main  engine  bells. 

B.         DIRECTIONAL  STABILITY  IN  THE  TWO  POINT  ATTITUDE 

Directional  stability  in  the  two  point  stance  is  different  than  the  directional  stability  of  the 
Orbiter  in  flight.  In  flight,  the  Obiter  would  typically  rotate  about  the  Center  Of  Gravity  (COG). 
In  the  two  point  stance,  assuming  no  lateral  skidding  of  the  mainwheel  tires,  the  Orbiter  is 
constrained  to  rotate  directionally  about  a  point  on  a  line  that  runs  laterally  through  the  areas  of 
contact  between  the  mainwheels  and  the  runway.  This  point  is  referred  to  as  the  Center  Of 
Rotation  (COR).  Directional  stability  is  governed  by  four  main  factors  in  the  two  point  attitude: 

(1)  The  relative  position  of  the  COG  to  the  COR. 

(2)  The  increased  friction  on  the  downwind  mainwheel  caused  by  the  rolling  moment 
generated  as  a  result  of  the  dihedral  effect  or  roll  due  to  sideslip,  Lp,  generated 
primarily  by  the  highly  swept  wing. 

(3)  The  aerodynamic  'weathercock'  effect  or  yaw  due  to  sideslip,  Np,  generated  primarily 
by  the  large  vertical  tail. 

(4)  The  relative  position  of  the  drag  chute  to  the  COR. 

The  first  factor  will  attempt  to  keep  the  Orbiter  traveling  in  a  straight  line  with  respect  to 
the  runway  in  a  zero  crosswind  condition  but  can  act  to  direct  the  Orbiter  away  from  a  crosswind 
as  a  result  of  the  second  factor.  The  third  factor  will  tend  to  fair  the  Orbiter  into  a  crosswind  in 
flight.  After  touchdown,  with  the  COR  shifting  aft,  this  'weathercock'  effect  becomes  a  much 
weaker  contributor.  The  forth  factor  will  attempt  to  fair  the  Orbiter  into  a  crosswind. 

1 .  Center  of  Gravity  (COG)  Effects 

The  nominal  COG  envelope  is  32  in.  long  and  3  in.  wide,  centered  laterally  and 
positioned  approximately  8  Vi  ft.  forward  of  the  mainwheels  [Ref  5].  The  COG  will  always  be 
well  ahead  of  the  mainwheels  in  the  two  point  stance.  The  COG  geometry  with  respect  to  the 
mainwheels  is  depicted  in  Figure  2. 


Nominal 

Center  of  Gravity       Mainwheels' 
Range 


Figure  2.  Orbiter  Center  of  Gravity  Location 


With  this  geometry,  the  COG  acts  as  a  pendulum  in  the  horizontal  plane.  The 
deceleration  of  the  Orbiter  caused  by  the  friction  in  the  mainwheels  and  the  main  wheel  brakes  acts 
on  the  COG  to  produce  a  force  in  the  Orbiter's  direction  of  motion.  This  force  is  analogous  to  the 
gravitational  force  on  a  normal,  vertically  oriented,  pendulum.  The  result  is  a  restoring  yaw 
moment  which  is  always  attempting  to  maintain  the  Orbiter  on  a  straight  track  in  the  direction 
from  the  COR  through  the  COG.  In  a  zero  crosswind  condition  with  no  mainwheel  braking,  the 
COR  will  most  likely  be  positioned  directly  between  the  mainwheels.  In  this  event,  the  Orbiter 
will  have  positive  directional  stability  in  the  direction  of  motion.  In  a  crosswind  case,  the  high 
dihedral  of  the  Orbiter  can  become  a  directionally  destabilizing  factor. 


2. 


Dihedral  Effect 


The  Orbiter  has  a  large  dihedral  effect  due  primarily  to  the  highly  swept  wing.  Dihedral  is 
the  tendency  for  an  aircraft  to  roll  away  from  sideslip.  In  a  landing  scenario,  the  sideslip  angle  can 
be  calculated  as  follows: 


/?=tan" 


crosswind    component 


Orb  iter     Indicated     Airspeed 

The  rolling  moment  generated  in  a  highly  swept  wing  aircraft  in  response  to  sideslip  can 
be  explained  by  two  factors.  First,  the  upwind  wing  senses  flow  closer  to  perpendicular  to  the 
leading  edge  which  decreases  flow  in  the  spanwise  direction;  thereby  increasing  lift.  The 
downwind  wing  has  the  opposite  effect  of  increased  spanwise  flow;  therefore  decreased  lift.  The 
second,  compounding  effect,  is  a  blanking  effect  on  the  downwind  wing  caused  by  the  large 
fuselage  structure  of  the  Orbiter.  This  reduction  in  airflow  also  decreases  lift. 

Another  contributor  to  positive  dihedral  is  the  high  vertical  tail.  The  center  of  pressure 
associated  with  the  airflow  across  the  vertical  tail  is  well  above  the  COR  which  would  cause  a 
further  roll  away  from  crosswind.  The  net  result  is  a  rolling  moment  away  from  crosswind  as 
depicted  in  Figure  3. 


^° 
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Left  Rolling 
Moment 
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"p^V^^       More 
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Side  Force  from 
High  Vertical  Tail 


Figure  3.  Orbiter  Dihedral 

This  rolling  moment  will  cause  an  increased  down  force  on  the  downwind  strut  and  a 
lighter  force  on  the  upwind  strut.  This  differential  load  will  cause  increased  friction  on  the 
downwind  wheel  and  a  subsequent  shift  in  the  directional  COR  to  the  downwind  side  of  Orbiter 
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centerline.   This  shift  in  the  COR.  in  conjunction  with  the  COG  effect  discussed  in  subsection  1, 
will  act  to  cause  the  Orbiter  to  deviate  away  from  crosswind  as  depicted  in  Figure  4. 


Figure  4.  Orbiter  Yaw  Due  to  Dihedral 


3. 


Weathercock  Effect 


The  weathercock  effect  is  typically  a  tradeoff  between  the  destabilizing  effect  of  the 
fuselage  and  the  stabilizing  effect  of  the  vertical  tail.  The  aft  location  of  the  COG  of  the  Orbiter 
in  flight  is  somewhat  destabilizing  since  the  stabilizing  moment  arm  associated  with  the  fin  is 
shorter  and  the  destabilizing  moment  arm  of  the  fuselage  is  longer.  When  the  Orbiter  attains  the 
two  point  stance,  the  COR  shifts  even  further  aft.  This  change  creates  a  condition  of  weak 
instability.  The  Orbiter  will  tend  to  point  away  from  a  crosswind  as  depicted  in  Figure  5. 


Left  Yawing 
Moment 


Vertical  Fin 
Force 


Figure  5.  Obiter  Yaw  Due  to  Weathercock 

4.         Drag  Chute  Effect 

The  drag  chute  was  incorporated  on  to  the  Orbiter  to  decrease  landing  rollout  distances  for 
certain  contingency  runways.  Along  with  the  decreased  rollout,  there  is  a  small  amount  of 
positive  weathercock  effect.  It  is  important  to  note  that  the  drag  chute  is  nominally  fully  disreefed 
three  seconds  prior  to  nosewheel  touchdown  as  described  in  section  II. A  [Ref.  2  and  3]. 
Therefore,  the  drag  chute  will  only  be  a  brief  factor  to  steady-state  directional  stability  in  the  two 
point  stance.  Yet  the  three  second  impulse  may  be  enough  to  perturb  the  steady  case.  The  drag 
chute  will  blow  downwind  from  a  crosswind.  Since  the  drag  chute  is  attached  aft  of  the  center  of 
rotation,  the  resulting  yaw  moment  on  the  Orbiter  is  into  the  crosswind  as  depicted  in  Figure  6. 
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Right  Yawing 
Moment 


Figure  6.  Obiter  Yaw  Due  to  Drag  Chute 

5.  Frequency  and  Damping 

A  summary  of  the  directional  stability  parameters  is  shown  in  Table  2  along  with  their 
relative  strengths. 

Table  2 .  Summary  of  Directional  Stability  Parameters 


Parameter 

Type 

Relative  Strength 

Direction 

COG 

Inertial 

Strong 

In  direction  of  motion 

Dihedral 

Aerodynamic 

Weak 

Away  from  crosswind 

Weathercock 

Aerodynamic 

Weak 

Away  from  crosswind 

Drag  Chute 

Aerodynamic 

Weak 

Into  crosswind 

All  of  these  factors  combine  to  yield  fairly  weak  directional  stability  in  the  two  point 
stance.  This  weak  stability,  when  modeled  with  a  simple  second  order  system,  would  yield  a  low 
open  loop  frequency  response  to  a  disturbance.  While  this  data  point  has  never  been  collected  in 
flight,  a  rough  estimate  of  what  one  might  call  'the  directional  natural  frequency',  o>d,  can  be  made 
by  observing  video  of  actual  closed  loop  cases  of  PIO.  The  videoThe  open  loop  natural 
frequency,  cod,  is  estimated  to  be  on  the  order  of  1  rad/sec.    The  only  yaw  damping,  Nr,  is  most 
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likely  due  to  the  large  vertical  tail  and  possibly  any  differential  friction  in  the  mainwheels.  Both 
of  these  factors  would  be  small  which  yields  very  low  directional  damping.  The  complete 
dynamic  stability  of  the  Orbiter  in  the  two  point  stance  can  be  modeled  by  the  following  second 
order  equation  where  NTd  is  total  damping  and  NTS  is  total  stability. 

s2  +  Nms  +  jVre 


where  coD  =  yjNj^ 


TTi 

and  the  damping  ratio,  C,D  =  - 

2o)D 

The  natural  frequency  and  damping  ratio  will  be  estimated  from  the  simulation  results.  A 
low  frequency  and  low  damping  ratio  will  indicate  the  tendency  of  the  FCS  to  be  susceptible  to 
PIO. 

C.  DIRECTIONAL  CONTROL  EN  THE  TWO  POINT  ATTITUDE 

Directional  control  in  the  two  point  stance  can  be  affected  with  two  controllers; 
aerodynamically  with  the  rudder,  and  mechanically  with  differential  braking.  Brakes  are  not 
employed  until  derotation  is  complete  to  avoid  excessive  nosewheel  slapdown  rates  and  to  reduce 
the  risk  of  a  blown  main  tire.  Braking  nominally  begins  at  140  knots  [Ref  3].  Differential 
braking  would  only  be  employed  in  extreme  situations,  such  as  a  flat  main  gear  tire.  Therefore, 
directional  control  in  the  two  point  stance  would  be  afforded  only  by  the  rudder. 

1.  Rudder 

The  rudder  is  actually  a  dual  purpose  control  surface  acting  as  both  a  directional  controller 
and  a  speedbrake  as  depicted  in  Figure  7. 
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31ft. 


Rudder  Neutral  -  Speedbrake  Closed 


Rudder  Right  -  Speedbrake  Open 


Figure  7.  Orbiter  Rudder/Speedbrake 

The  surface  is  mechanically  limited  to  a  maximum  deflection  rate  of  16  deg/sec  normally 
[Ref.  6].  Postlanding  the  deflection  rate  is  further  limited  tol3deg/sec  since  the  speedbrake  is  full 
open.  The  software  limits  the  deflection  rate  to  14  deg/sec  so  that  sufficient  hydraulic  fluid  is 
available  to  other  control  surfaces.  In  the  event  of  the  lost  of  one  or  two  hydraulic  system,  the  rate 
is  limited  to  12  deg/sec  and  7  deg/sec,  respectively  [Ref.  4].  The  surface  areas  is  approximately 
97  ft.2  and  can  generate  up  to  2.25  lateral  g's. 

2.  Directional  FCS  Susceptibility  to  PIO 

The  FCS  has  many  limitations  which  result  in  the  high  potential  for  directional  PIO  in  the 
two  point  stance.  The  slow  surface  rate  of  the  rudder,  slow  response  time,  high  control  power  of 
the  rudder  due  to  the  large  surface  area,  and  rudder  placement  with  respect  to  the  gear  are 
conducive  to  PIO.  If  a  large  pedal  input  is  sustained  until  the  yaw  rate  is  well  underway,  loss  of 
control  is  probable.  The  slow  surface  rate  does  not  allow  an  opposite  input  to  be  effected  quicky 
enough  to  counteract  the  yaw  rate. 

In  addition,  ARC  VMS  studies  have  indicated  rudder  rate  saturation  is  encountered  in 
high  crosswind  landings  or  in  performing  postlanding  lateral  maneuvers  stress  tasks  [Ref.  7]. 
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III.  RUDDER  CONTROL 

This  chapter  describes  the  rudder  control  path  including  the  software  command  path  and 
rudder  actuator.  Detailed  drawings  and  description  of  the  Orbiter  FCS  control  laws  are  found  in 
the  Software  Requirements  Document  [Ref  8]  and  Lockheed  drawings  [Ref  9]. 

A.         RUDDER  COMMAND  PATH 


The  control  logic  for  the  rudder  surface  may  be  broken  down  into  two  basic  steps.  The 
first  converts  rudder  pedal  deflection  angle  into  lateral  acceleration  error  command.  The  Orbiter's 
flight  control  system  uses  lateral  acceleration  feedback  to  null  errors  to  the  commanded  lateral 
acceleration  and  provide  directional  stability  augmentation.  The  second  step  converts  the  lateral 
acceleration  error  command  into  a  yaw  rate  which  is  filtered  of  Orbiter  structural  bending  modes 
and  finally  converted  into  a  command  to  the  control  surface  for  degrees  of  deflection.  This 
command  is  sent  to  Priority  Rate  Limiting  (PRL),  which  limits  the  rudder  deflection  rate 
depending  on  the  number  of  hydraulic  systems  available.  These  two  steps  will  be  henceforth 
referred  to  as  the  rudder  command  feedback  and  rudder  deflection  command  portions. 


1. 


Rudder  Command  Feedback 


A  simplified  block  diagram  of  rudder  pedal  position  conversion  to  lateral  acceleration 
error  command  is  shown  in  Figure  8. 
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Figure  8.  Rudder  Command  Feedback 
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The  various  components  of  the  rudder  command  feedback  path  are  discussed  in  the 
subsections  below. 


a. 


Rudder  Pedals 


There  is  a  Rudder  Pedal  Transducer  Assembly  (RPTA)  located  at  both  the 
Commander  and  Pilot  stations.  Each  RPTA  consists  of  two  pedals  and  a  transducer  which 
converts  degrees  of  rudder  pedal  deflection  into  electronic  signal.  The  RPTA  uses  a  time  constant 
of  100  ms  in  the  mechanical  to  electronic  transfer.  This  signal  is  processed  at  6.25  Hz  and  passed 
to  the  rudder  command  shaper  with  the  name  tag  DRMAN. 

b.         Rudder  Command  Shaper 

The  rudder  command  shaper  converts  DRMAN  into  a  shaped  and  scaled  rudder 
pedal  command  named  DRMS  as  shown  in  Figure  9. 


DRMAN 


22  5  max 


22.5  max 


DRMS 


Figure  9.  Rudder  Command  Shaper 

The  shaping  schedule  is  parabolic  and  conforms  to  the  following  function: 
DRMS  =  (0.042)DRMAN2  +  (0.\3\)DRMAN .  A  maximum  deflection  limit  of  22.5  degrees 
is  imposed  on  the  output  signal. 

c.  Rudder  Pedal  Gain  Scheduler 

The  rudder  pedal  gain  scheduler  produces  a  gain  signal,  which  when  multiplied 
with  the  shaped  and  scaled  rudder  pedal  command,  DRMS,  produces  a  lateral  acceleration 
command.  This  gain  is  scheduled  as  a  function  of  True  Airspeed  (TAS)  and  is  given  the  name 
GNYDRM  as  depicted  in  Figure  10.  This  gain  also  converts  the  units  of  the  signal  from  degrees 
of  pedal  deflection  to  lateral  g's.  The  rudder  pedal  gain  scheduler  samples  at  1 .04  Hz. 
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Figure  10.  Rudder  Pedal  Gain  Scheduler 


d. 


Rudder  Trim  Command 


The  rudder  trim  command  signal  is  generated  by  converting  the  command  signal 
from  the  rudder  trim  panel  into  a  commanded  lateral  acceleration  bias.  This  is  performed  with  a 
rectangular  integrator  which  ramps  the  output  signal  and  also  limits  the  maximum  bias  to  0. 1  g's. 
This  signal  is  sampled  at  6.25  Hz.  The  output  signal  is  the  lateral  acceleration  trim  signal,  DRTI, 
which  is  summed  with  GNYDRM  from  the  rudder  pedals  to  produce  a  total  lateral  acceleration 
command,  DRTMS,  as  depicted  in  Figure  1 1 . 
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Figure  1 1 .  Rudder  Trim  Commands 
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e. 


Lateral  Accelerometers 


There  are  four  Accelerometer  Assemblies  (AA)  located  forward  of  the  crew 

compartment.     These  AA's  directly  measure  lateral  acceleration  and  simultaneously  filter  the 

signal  as  shown  below: 

Measured     Acceleration 

r-  — »  Lateral     Acceleration    Signal(NYP) 


15' 


( 


146' 


s2  +2(.5)(\5)s  +  \52As2  +2(.8)(l46)s  +  146: 


compensator. 


This  signal  is  sampled  at  25  Hz  and  passed  to  the  lateral  acceleration  feedback 


/  Lateral  Acceleration  Feedback  Compensator 


The  lateral  acceleration  feedback  compensator  compares  the  actual  lateral 
acceleration,  NYP,  as  passed  from  the  AA's,  to  the  commanded  lateral  acceleration,  DRTMS,  to 
generate  an  error  signal,  or  additional  lateral  acceleration  command.  A  low  pass  filter  is 
incorporated  in  this  procedure  to  smooth  the  output.  The  resulting  signal  is  sampled  at  25  Hz. 
The  process  is  depicted  in  Figure  12. 
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Figure  12.  Lateral  Acceleration  Feedback  Compensator 


2. 


Rudder  Deflection  Command 


A  simplified  block  diagram  of  lateral   acceleration   conversion  to  rudder  deflection 
command  is  shown  in  Figure  13. 
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Figure  13.  Rudder  Deflection  Command 


(3    Multiplier 
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The  various  components  of  the  rudder  deflection  command  path  are  discussed  in  the 
subsections  below. 


a. 


Lateral  Acceleration  Gain 


The  lateral  acceleration  gain  scheduler  produces  a  gain  signal,  GRAY,  which 
when  multiplied  with  the  lateral  acceleration  error  command,  DRTMS,  produces  a  yaw  rate 
command,  DAY.  This  gain  is  scheduled  as  a  function  of  Mach  as  depicted  in  Figure  14.  This 
gain  also  converts  the  units  of  the  signal  from  g's  to  deg/sec.  The  sample  rate  for  this  scheduler  is 
1.04  Hz. 
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Figure  14.  Lateral  Acceleration  Gain  Schedule 


Cosine  of  the  Angle  of  Attack  (a) 


A  dimensionless  signal,  cos(a),  is  mixed  with  the  yaw  rate  command  signal  and 
summed  with  the  yaw  rate  error  to  produce  a  yaw  rate  command  about  an  inerrial  axis.  This 
signal  will  always  reduce  the  yaw  rate  command  signal  in  producing  a  yaw  rate  with  respect  to  the 
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runway.  The  cos(a)  is  sampled  at  1 .04  Hz.   The  output  signal  from  this  mix  is  referred  to  as  the 
unfiltered  yaw  rate  command,  DRRCUF. 

c.  Rate  Gyro  Assemblies 

There  are  four  Rate  Gyro  Assemblies  (RGA's)  located  on  the  aft  bulkhead  below 
the  floor  of  the  payload  bay.  Each  of  the  RGA's  sense  rotation  about  the  three  vehicle  axes  and 
simultaneously  filter  the  signal  as  shown  below: 

Measured     Acceleration 


50 


s  +50 


-—, ^^Body    Rates    (P,R,Q) 


control. 


V +2(.7)(420Xv  +  4202J 

These  rates  are  the  primary  feedback  to  the  flight  control  system  to  maintain 

d.  Yaw  Rate  Structural  Filter 

A  second-order  structural  bending  filter  is  applied  to  the  yaw  rate  feedback. 

e.  Yaw  Rate  Error  Determination 

The  filtered  body  yaw  rate,  DRPRM,  is  summed  with  the  product  of  the  sideslip 
angle  times  the  cosine  of  the  angle  of  attack,  BCOS  ALF,  and  yaw  rate  command,  DAXFDC,  to 
determine  the  yaw  rate  error,  DRERR,  as  shown  in  Figure  15.  The  sideslip  angle  is  not  used  in 
the  approach  and  landing  flight  phase.  Therefore,  BCOSALF  is  zero  in  this  analysis.  The 
resulting  signal  is  sampled  at  25  Hz. 

BCOSALF     x-fci      DAXFDC 


DRPRM 


DRERR 
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Figure  15.  Yaw  Rate  Error  Determination 
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/  Structural  Bending  Filter 

A  structural  bending  filter  is  included  in  series  to  remove  any  known  Obiter  body 
bending  modes.  The  resulting  signal  is  the  yaw  command,  DRRC,  which  is  processed  at  25  Hz. 
The  filter  is  made  of  two  consecutive  second  order  notch  filters.  This  fourth  order  system  will 
inherently  result  in  transport  delays  which  would  only  contribute  to  the  time  lag  problem. 


g- 


Yaw  Rate  Gain 


The  yaw  rate  gain  scheduler  produces  a  gain  signal,  GDRC.  This  gain  and  a 
constant,  GERYFWDS,  is  multiplied  with  the  yaw  command  signal,  DRRC,  to  produce  an  initial 
rudder  deflection  command,  DRCPF.  This  gain  is  scheduled  as  a  function  of  Mach  and  dynamic 
pressure  as  depicted  in  Figure  16.  This  gain  also  converts  the  units  of  the  signal  from  deg/sec  to 
deg.  The  sample  rate  for  this  scheduler  is  1 .04  Hz. 


Dynamic 

Pressure 

(psf) 

Figure  16.  Yaw  Rate  Gain  Scheduler 

h.  Rudder  Deflection  Limiter 

The  rudder  deflection  limiter  provides  a  maximum  rudder  deflection  limit  in  the 
positive  and  negative  directions.  The  deflection  limit,  DRC  LIM,  is  24. 1°  and  27. 1°  in  flight  and 
post  Weight  On  Wheels  Lock-On  (WOWLON),  respectively.  Within  these  limits,  the  limiter 
passes  a  one-to-one  correlation  from  the  filtered  rudder  deflection  command,  DRCPF,  to  the 
actual  rudder  deflection  command,  DRC.  The  WOWLON  signal  is  sampled  at  6.25  Hz,  but  the 
rudder  deflection  limiter  processes  at  25  Hz.  The  system  is  configured  as  shown  in  Figure  1 7. 
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Figure  17.  Rudder  Deflection  Limiter 


/.  Priority  Rate  Limiting  (PRL) 

The  rudder  deflection  command  is  sent  to  PRL.  PRL  limits  control  surface 
deflection  rates  to  manage  the  amount  of  hydraulic  fluid  sent  to  each  of  the  actuators.  There  are 
three  Auxiliary  Power  Unit  (APU)/hydraulic  systems  available  to  provide  hydraulic  fluid  to  the 
control  surface  actuators.  Normally,  with  three  hydraulic  systems  available,  PRL  limits  the  rudder 
deflection  rate  to  14  deg/sec  [Ref  4].  This  slow  deflection  rate  is  conducive  to  the  pilot's  rudder 
inputs  being  out  of  phase  with  Orbiter  response.  The  deflection  rate  is  further  reduced  depending 
upon  the  number  of  hydraulic  systems  failed.  With  one  and  two  systems  failed,  the  rudder 
deflection  rate  is  limited  to  1 2  and  7  deg/sec,  respectively. 


B. 


RUDDER  ACTUATOR 


The  rudder/speedbrake  actuators  receive  hydraulic  fluid  from  all  three  hydraulic  systems 
through  a  switching  valve.  The  switching  valve  produces  one  hydraulic  output  pressure  source  for 
the  rudder/speedbrake  servovalves.  The  power  drive  unit,  mechanical  actuator,  and  actuator  drive 
shaft  are  shown  in  Figure  18. 
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Figure  18.  Rudder/Speedbrake  Actuator  Location  [Ref.  6] 
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IV.  EQUATIONS  OF  MOTION  AND  STABILITY  AND  CONTROL 

DERIVATIVES 

A.         ASSUMPTIONS 

The  following  assumptions  were  used  in  developing  the  Equations  of  Motion  (EOM)  for 
the  Obiter  in  the  two  point  stance: 

1 .  Inertial  frame  fixed  to  the  flat  Earth  since  the  Orbiter  is  on  the  ground. 

2.  Constant  vehicle  mass  properties  since  there  is  no  propellant  consumption  and 
nothing  ejected  until  drag  chute  deploy. 

3.  Vehicle  is  symmetric  about  the  XZ  plane.  Therefore,  Ixy  =  Iyz  =  Cyo  =  C\0  =  Cno  =  0. 

4.  Equations  of  motion  are  linear  about  a  steady  state  flight  condition.  Small 
perturbation  theory  implies  that  the  squares  and  products  of  perturbations  may  be 
neglected.  For  every  perturbation  angle  a,  Sin  (a)  =  a  and  Cos  (a)  =  1 . 

5.  Atmospheric  property  variations  are  negligible.  Air  density,  temperature,  pressure 
changes  are  negligible  over  a  twenty  second  period  and  less  than  10,000  feet  rolling 
distance  on  the  ground. 

6.  Steady  state  rotation  rates,  side  velocity,  and  bank  angle  are  all  zero.  The  Orbiter  in 
the  two  point  stance  will  not  perform  rapid  maneuvers  or  fly  with  significant  steady 
state  sideslip. 

7.  The  air  flow  is  quasi-steady.  In  other  words,  the  pressure  distributions  adjust  instantly 
to  rate  changes.  This  assumption  is  generally  valid  for  angle  of  attack  rate  a  if  the 
following  condition  is  satisfied  [Ref.  10]:  ac/2Ui  <  0.04.  For  the  Orbiter,  the  chord 
length  c  =  39.56  ft.  and  the  initial  velocity  Ui  =  353  ft./sec.  The  inequality  is 
satisfied  if  a  <  0.714  rad/s  =  40.9  deg/s.  In  the  two  point  stance  prior  to  the  initiation 
of  derotation,  a  is  essential  constant.  After  initiation  of  derotation,  a  is  1.5  deg/sec. 
The  stability  derivatives  with  respect  to  a  may  be  neglected. 

8.  Morion  in  the  Z  direction  is  constrained  due  to  the  runway.  Therefore,  the  downward 
velocity  W=0  and  the  associate  forces  and  moments  are  zero. 

9.  Vehicle  is  at  constant  a.  Therefore,  all  pitching  moment  derivatives  are  zero,  the 
pitch  rate,  q,  is  zero,  and  all  control  derivatives  with  respect  to  q  are  zero. 
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10.  The  elevons  are  at  trail  (0°)  at  touchdown  and  remain  there  until  the  initiation  of 
derotation  [Ref.  2].  Therefore,  all  variations  of  the  coefficients  due  to  elevator  and 
aileron  are  zero. 

1 1 .  Forward  velocity,  U  is  constant. 

B.         EQUATIONS  OF  MOTION 

The  general  EOM  in  the  body- fixed  coordinate  system  obtained  from  Roskam  [Ref.  10] 
are  defined  below.  The  subscript  denotes  the  axis  and  the  apostrophe  denotes  time  derivative.  All 
symbols  are  defined  in  the  List  of  Symbols,  Acronyms,  and/or  Abbreviations  on  page  xii. 
Although  the  control  problem  under  investigation  is  a  directional  control  problem,  both  the 
longitudinal  and  lateral/directional  equations  are  included  in  the  derivation  and  simulation  for 
completeness.  The  general  force,  moment,  and  kinematic  equations  are: 
Fx  =  m(U  -VR+WQ)+mgSin(0) 
Fy  =  m(V  +UR-WP)-mgCos(0)Sin(O) 
F2  =  m(W-UQ+VP)-mgCos(O)Cos(0) 

L  =  IxxP'-I«R'-lKPQ+(Izz-Iyy)RQ 

M  =  IyyQ  +(IXX-IZZ)PR+IXZ(P2-R2) 
N  =  InR  -IJP  +(Iyy-Ixx)PQ+IxzQR 
P  =  O-v|/Sin(0) 
Q  =  0  Cos(0)+i|/  Cos(0)Sin(4>) 
R  =  V|/  Cos(0)Cos(<D)-0  Sin(O) 

Each  variable  is  the  sum  of  a  steady  state  components,  denoted  with  a  subscript  "  1 ",  and 
perturbation  components,  denoted  with  lower  case: 
U  =  U,+u  V  =  V,+v  W=W,+w 

P  =  P,+p  Q  =  Qi+q  R  =  R,+r 

4'  =  4'1+ij/  0  =  0,+6  0  =  <D1+<J) 

T  =  ^,  +vj/  0  =  0,  +0  <P  =  O,  +((> 
Fx+Fxl  +  fx  Fy  +  Fyl  +  fy  Fz+Fzl  +  fz 
L  =  L,+1  M  =  M,+ma        N=N,+n 

The  steady  state  forces,  moment,  and  kinematic  equations  are  as  follows: 
Fxl  =  m(U  rV.R.+W.QO+mgSinC©,) 
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Fy]  =  m(V'i+UiRrWiPi)-mgCos(0i)Siii(a>i) 
Fzl  =  m(W  rU^^V.PO-mgCosCOOCosiG,) 

L,  =  I.xP  rlxzR  l-IxzPlQ.+(Izz-Iyy)RlQl 

M,  =  lyyQ ,  +(Ixx-IZz)P.Ri+UP.2-R12) 
N,  =  I^R ,  -U>',  +(I>y-Ixx)P.Q.  +1xzQiRi 
Pi  =  O'i  -\\i  ,Sin(0) 
Q,  =  0',Cos(<Di)+v|/  ^Cos^Sin^,) 
R,  =  h;  ,Cos(0!)Cos(Oi)-0 ^SinCOO 

The  pitch  and  bank  angles  are  defined  such  that: 
Sin(9)  =  0  Sin((|>)  =  <j>  Cos(9)  =  Cos(<|>)  =1.0 

This  small  angle  assumption  restricts  the  pitch  and  bank  angle  perturbations  to  less  than 
15  degrees.    This  limit  is  reasonable  for  landing  and  rollout  and  does  not  constitute  any  serious 
restriction.     Perturbations  are  assumed  to  be  sufficiently  small  such  that  products  and  cross- 
products  of  the  perturbations  are  negligible  with  respect  to  the  perturbation  itself: 
pv  =  pw  =  pq  =  pr  =  qu  =  qw  =  qr  =  ru  =  rv  =  0 
u2  =  v2  =  w2  =  p2:=r2  =  (p9  =  i|/9  =  vj/(p  =  (p9  =0 

The  perturbation  force,  moment,  and  kinematic  equations  are  as  follows: 
fx  =  m(g8Cos(0)+Wiq-ViH-u'-RiV+Qiw) 

fy  =  mC-g^CosCOOCosCOO+gGSinCOOSinCOO-W^+U.r+Riu+v  -P,w) 
fz  =  mCgipCosCOOSinCOO+geCosCOOSinCO^Vip-Uiq-Qiu+Piv+w") 
1  =  IxxP  -(QiP+Piq)Ixz+(-Iyy+Izz)(Riq+Q.r)-lxzr 
ma  =  2IxzPip+(Ixx-Izz)R1p+Iyyq  +(Ixx-Izz)Pir-2IxzR,r 
n  =  (-Ixx+Iyy)(QiP+Piq)-IxzP  +(Riq+Qir)IxZ+IZ2r 
p  =  -94*  ,Cos(0,)  -  i|/'Sin(0i)+q>" 

q  =  9  CosCOO+H/CosCOOSmCOO^T  ,Sin((I)1)Sm(01)+4/ jCosC^OCosC©,)-© ^SmCO,)^ 
r  = -(0  ^CosCOO+T^CosCOOSmCO,))^^  Cos(OOCos(0,)-9Sin(O,)-9T  ^CosCOOSmC©,) 

1.  Longitudinal  Equations  of  Motion 

The  longitudinal  forces  in  the  stability  axis  system  are  defined  as: 
Fxl  =  -D     Fzi  =  -L     D,  =  CD,qdynS     Lj  =  CLiqdynS     M,  =  CmiqdynSc 

The  non-dimensional  stability  and  control  derivatives  in  the  stability  axis  are  defined  as: 
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d  d  d  d 

c    = c      c    = — c      c c      c     = — c 

<~DU  u    ^D         ^Da  fa  ^D         ^Dq    -  qc    ^D         ^DSt     "   gg    ^D 


Ux  W 


L»  "  a  u  Cl      Cl°  ~  daCl  Cl<  ~  sqcCl      Cls>  ~  d8e     L 
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d                       _^_  d                           d 
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The  longitudinal  dimensional  stability  axis  derivatives  are  defined  as: 

Axial  force: 

x  _-^S{Cdu+2Cd)    1  _~<l^s{cDa+2CL)       ft                 ^x^   i 


mUl  sec  m  rad  sec  Ux  sec 

-q+frCp,    ft  -q<*y»scnSt     ft 

A„   —  _        __  .  Aj     - 


lq  2mUx       rad  sec  s'  m         rad  sec2 

Normal  force: 


Z„    1 


mf/j  sec  /w  rad  sec2  Ux  sec 

z  _-q+&c^    ft       z  _~^scLSe     ft 

q  2mUx       rad  sec         ^  w  rad  sec2 

Pitching  moment: 

m..  = — — ; M„  = : t      M...  = 


1„UX  ft  sec       "a  I„       sec2  Ux  //sec 

q^c2c     1  qdynScCm     i 


2/w,t/,      sec  *  /„        sec2 

Recall  that  motion  in  the  Z  direction  is  constrained  due  to  the  runway.  Therefore,  the 
downward  velocity  W=0  and  the  associate  forces  and  moments  Xw,  Zw,  and  Mw  are  zero  and  not 
included  in  the  state  space  equation. 

The  pitch  moment  of  inertia  changes  slightly  when  the  Orbiter  touches  down.  The  center 
of  rotation  shifts  from  the  center  of  mass  to  the  axis  formed  by  the  line  between  the  points  where 
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the  main  gear  tires  meet  the  runway.  The  parallel  axis  theorem  [Ref.  1 1]  allows  a  simple  shift  of 


the  moment  of  inertia  with  the  following  equation:    /       =  /     +  mh"  where  lymg  is  the  moment  of 

inertia  about  the  main  gear,  m  is  the  mass  of  the  Orbiter,  and  h  is  the  distance  between  the  center 
of  mass  and  the  new  rotation  point.  All  pitch  moment  dimensional  stability  derivative  equations 
will  be  adjusted  accordingly  by  multiplying  by  Iyy/Iymg. 

The  following  are  the  perturbed  longitudinal  equations  of  motion: 

u'=  Xuu  +  Xaa  +  Xqq  -  gcos(0>+  XsSe 
,     Zmu     Zaa     {Ut+Z^     gsin(©>     ZsSe 
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2.  Lateral-Directional  Equations  of  Motion 

The  lateral-directional  forces  are  defined  as: 
Fyl  =  CyqdynS     L)  =  Cnqd>T)Sb     Nj  =  CnqdynSb 

The  non-dimensional  stability  and  control  derivatives  are  defined  as: 
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d  d  d  d  d 

c   = — c     c   = c     c  = c     c    = — C      C    - — c 

ne     ap*      ""     ^pp  r       it  -     d8a  r     d8r 

The  lateral-directional  dimensional  stability  and  control  derivatives  are  defined  as: 
Side  force: 

Y  =  q+nSCy,     ft        y  =  q^sbcyp    ft       ^  =  i+&e,f    ft 


m       rad  sec2  2mUx     rad  see  2mU{     rad  see 


<J+.SCys.  ft  J^ys,  ft 

I  x     —  -i  Is     — 


&a  m        rad  sec2        Sr  m        rad  see' 

Rolling  moment: 


q+,sbch    i  g+,sb2clp  l         _q^sb2clr  1 

J-*  ft    —  _  -7  A-/  „    —  _   _       __  /-/_    — 


2 


sec'  2/^7,     sec        r        21xxUl     sec 


qdynSbc,s     x  q^sbc, 


I„        sec2  L*  !„        sec2 

Yawing  moment: 

Ar     q^sbcnp    i  g^c     i      a,     ^^2c„r   i 

N«  = ; r-  TV    = — — - — N. 


Ia        sec2  p  '      2/„I/,      sec         r         2/sC7 


sec 


s-  Ia         sec2  *-  /„         sec2 

The  yaw  moment  of  inertia  changes  slightly  when  the  Orbiter  touches  down.  The  center 
of  rotation  shifts  from  the  center  of  mass  to  a  point  between  the  main  gear.  It  is  assumed  that  this 
point  will  be  directly  between  the  main  gear.   The  parallel  axis  theorem  [Ref  1 1  ]  allows  a  simple 

shift  of  the  moment  of  inertia  with  the  following  equation:    /      =  I „  +  mh  where  Izmg  is  the 

moment  of  inertia  about  the  main  gear,  m  is  the  mass  of  the  Orbiter,  and  h  is  the  distance  between 
the  center  of  mass  and  the  new  rotation  point.  All  yaw  moment  dimensional  stability  derivative 
equations  will  be  adjusted  accordingly  by  multiplying  by  I^/I^g- 

The  roll  moment  of  inertia  is  assumed  to  remain  referenced  to  the  center  of  mass  as  the 
main  gear  struts  are  free  to  compress  and  extend  which  permits  the  Orbiter  to  roll  about  a  point 
near  the  center  of  mass.  1^  will  not  be  adjusted. 
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The  following  are  the  perturbed  lateral-directional  equations  of  motion: 
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For  small  0's  and  small  1^,  inertial  coupling  is  negligible;  therefore'  the  I^r'/I^  and  I^pVI^  terms 
are  disregarded. 

In  state  space  form,  X,=AX+BU: 
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3.    Longitudinal  and  Lateral-Directional  EOM  in  State  Space  Form 

The  full  8x8  state  space  equation  is  shown  below:  X'=AX+BU.  The  longitudinal  and 
lateral-directional  EOM  are  decoupled  in  this  analysis.  The  longitudinal  EOM  need  not  be  solved 
to  analyze  this  problem.  The  longitudinal  EOM  are  included  for  completeness. 
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4.  Determination  of  the  Non-Dimensional  Stability  Derivatives 

The  non-dimensional  force  and  moment  coefficients  were  determined  in  wind  tunnel 
testing  and  updated  based  on  actual  flight  data.  The  equations  and  coefficients  are  found  in  the 
Operational  Aerodynamic  Data  Book  Volumes  I  and  III  [Ref  12]. 

a.  Longitudinal  Stability  Derivatives 

The  Taylor  series  first  order  expansion  of  aerodynamic  coefficients  is  given  by: 

Cdi  =  Cdo  +  CdoOI  +  Cd&{&)\ 
Cl\  =  Clo  +  Clo£1\  +  Cl&(  &)\ 
C./wi  —  Crao  +  0-ma£7l  +  Cm<&\Ck)\ 

The  elevons  are  at  trail  (0°)  at  touchdown  [Ref.  2].  Postlanding,  the  elevon  is 
deflecting  as  required  to  balance  the  longitudinal  moment  equation  and  maintain  a  constant  pitch 
attitude.  However,  elevon  deflection  is  minimal  until  the  start  of  denotation.  At  this  time,  the 
elevons  ramp  quickly  downward  as  theta  decreases  at  a  rate  of  1-2  deg/sec.  Therefore,  (5e)i  is 
negligible  in  the  two  point  stance  and  set  to  0  in  these  three  equations. 

The  coefficients  [Ref.  1 2]  are  defined  as: 

M)l     :"    CqbaSIC    +AC()high    alt    "'"ACoreai    gas    +ACoviscous    mteraction    +AC()BF    +AC()elevator/aileron    +  AC( )p|unging 
+AC()pltch  rate  +AC()LG  +AQ)GE  +AC()sjits  pod  +AC(  Vehicle 

where  the  subscript  ( )  denotes  lift  force  (L),  drag  force  (D),  or  pitching  moment  (m)  and: 

Cqbasic  =  Basic,  full-scale,  freestream,  rigid  body  force  or  moment  coefficient. 

AC()h,ghait=  Change  in  force  or  moment  coefficient  due  to  high  alt.  effects  (h>285,000  ft). 

AC()reai  gaS  =  Change  in  force  or  moment  coefficient  due  to  real  gas  effects. 

ACoviscous  interaction =  Change  in  force  or  moment  coefficient  due  to  viscous  interaction  (VK  effects). 

AC()bi  =  Change  in  force  or  moment  coefficient  due  to  body  flap  deflection. 

AC( ^levator  aileron  =  Change  in  force  or  moment  coefficient  due  to  elevator  and  aileron  deflection. 

AC0piunging  =  Change  in  force  or  moment  coefficient  due  to  rate  of  change  of  alpha. 

ACopltCh  ratc  =  Change  in  force  or  moment  coefficient  due  to  pitch  rate. 

AC()LG  =  Change  in  force  or  moment  coefficient  due  to  fully  extended  landing  gear. 

AC()ge=  Change  in  force  or  moment  coefficient  due  to  proximity  of  the  ground. 

AC()Sjitspod=  Change  in  force  or  moment  coefficient  due  to  addition  of  silt  pods. 
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AC(  (vehicle =  Change  in  force  or  moment  coefficient  due  to  geometry  differences  between  vehicles. 

AC()h,ghait  and  AC()rea|gas  are  not  applicable  below  an  altitude  of  265,000  ft.  and  M 
<  6.0,  respectively.  AC0viscous  interaction  is  not  a  factor  since  V*  is  zero.  AC()e|evator,aiieiun  is  assumed  to 
be  zero  since  elevator/aileron  deflection  is  considered  negligible.  AC()piunging  and  AC()pltch  rate  are 
not  a  factor  with  the  Orbiter  on  the  ground  in  the  two  point  stance  since  the  pitch  rate  and  angle  of 
attack  are  essential  constant  until  the  initiation  of  derotation.  This  analysis  will  assume  Orbiter 
Vehicle  (OV)  103  which  does  not  have  a  silts  pod;  therefore'  AC()slits  pod  is  not  applicable.  The 
above  equation  can  be  reduced  to: 

C()I  =  Q)BASlC+AC()BF+AC()LG+AC()GE++AC()vehIC|e 

Each  of  the  terms  in  the  above  equation  was  obtained  by  linear  interpolation  of 
data  extracted  from  a  table  at  a  =  10°,  M=0.4  and  a  =  10°,  M=0.25  to  obtain  data  applicable  for 
M=0.325  (353  fps)  for  both  the  lift  and  drag  coeffiecients.  da  was  determined  by  breaking  each 
of  Cqbasic  AQ)bf,  AC()lg,  AC()Ge,  and  AC()vehiCie  terms  into  do  and  CLa  components  and  summing 
the  Clo  contributions.  Cr**  was  determined  in  the  same  manner. 

Since  U  is  considered  constant  in  this  analysis,    CD  and  and    CL  are  zero. 

CD  and  CL    are  zero  since  the  pitch  rate  is  zero.   The  lift  and  drag  coefficients  pertinent  to  this 
analysis  are  summarize  in  Table  3. 

Table  3.  Lift  and  Drag  Coefficients 


Lift  and  Drag  Coefficients 

Value  (per  °) 

Cl\ 

0.54055 

Cux 

0.05948 

Cd, 

0.16712 

Cdo 

0.00992 

b.  Lateral-Directional  Stability  Derivatives 

The  Taylor  series  first  order  expansion  of  the  lateral-directional  aerodynamic 
coefficients  are  provided  for  completeness.  The  Taylor  series  first  order  expansion  of 
aerodynamic  coefficients  is  given  by: 

Cy\  =  Cyo  +  Cyp/3\  +  Cyda\8a)\  +  C  y<>{&)] 
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Ch  =  Clo  +  Cipfi\  +  Cl6a{da)\  +  ClSr{5r)\ 
Cm  —  Cno  +  Cnfifll  +  Cn&\da)\  +  Cn&\5r)\ 

The  lateral-directional  bias  coefficients  are  assumed  to  be  zero  for  a  symmetrical 
airframe:  Cyo  =  Ch  =  Cno=  0.  At  touchdown,  a  software  flag  denoted  'flat  turn1  is  set  [Ref  13]. 
When  this  flag  is  set,  lateral  motion  is  accomplished  without  elevon  motion,  using  rudder  only. 
Therefore,  aileron  deflection  is  negligible  in  the  two  point  stance  and  (&)i  are  set  to  zero  in  the 
above  three  equations.  The  coefficients  [Ref  1 2]  are  defined  as: 

M)l  =  M)sideslip  +AC()alieron  +AC()rlldder  +AC()rol]  rate  +AC()yaw  rate  +AC()sjlts  pod  +AC()vehicle 

where  the  subscript  ( )  denotes  side  force  (y),  yawing  moment  (n),  or  rolling  moment  (1),  and: 
Q)sidesiiP=  Basic,  full-scale,  freestream  force  or  moment  coefficient  due  to  sideslip  angle. 
AC()ajieron =  Change  in  force  or  moment  coefficient  due  to  aileron  deflection. 
AC()rudder=  Change  in  force  or  moment  coefficient  due  to  rudder  deflection. 
AC()roii  rate  =  Change  in  force  or  moment  coefficient  due  to  roll  rate. 
AC()yaw  rate =  Change  in  force  or  moment  coefficient  due  to  yaw  rate. 
AC()sl|ts  pod  =  Change  in  force  or  moment  coefficient  due  to  addition  of  silt  pods. 
AC()vehicie  =  Change  in  force  or  moment  coefficient  due  to  geometry  differences  between  vehicles. 
AC()ajieron  may  be  dropped  from  this  equation  since  aileron  deflection  is  negligible. 

Each  of  the  lateral-directional  stability  and  control  derivatives  was  obtained  by 
linear  interpolation  of  data  extracted  from  a  table  at  a  =  10°,  M=0.4  and  a  =  10°,  M=0.25  to 
obtain  data  applicable  for  M=0.325  (353  fps).  The  side  force,  rolling  moment,  and  yawing 
moment  coefficients  pertinent  to  this  analysis  are  summarize  in  Table  4. 
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Table  4.  Side  Force,  Rolling  Moment,  and  Yawing  Moment  Coefficients 


Side  Force,  Rolling  Moment, 

and  Yawing  Moment 

Coefficients 

Value  (per  °) 

c» 

-0.01705 

c,. 

0 

Cy, 

0 

C,, 

0.00138 

(\ 

-0.00205 

C\ 

-0.3189 

Q 

0.1890 

C<„ 

0.0032 

C", 

0.0010 

C», 

0.0280 

Q 

-0.2925 

c«. 

-0.00054 
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V.  PROPOSED  FCS  MODIFICATIONS  TO  IMPROVE 
DIRECTIONAL  CONTROL  IN  THE  TWO  POINT  STANCE 

The  rudder  FCS  susceptibility  to  PIO  results  from  the  slow  surface  rate,  slow  response 
time  from  rudder  pedal  input  to  rudder  deflection,  the  power  of  the  rudder,  and  the  placement  of 
the  rudder  with  respect  to  the  gear.  Increasing  the  surface  rate  requires  hardware  modification. 
Analysis  of  hardware  modifications  are  beyond  the  scope  of  this  thesis.  Transport  delay  might  be 
decreased  by  shortening  time  delays  in  the  FCS,  linearizing  the  rudder  pedal  command  shaper, 
and/or  disabling  Ny  feedback.  The  power  of  the  rudder  might  be  reduced  by  decreasing  the 
rudder  pedal  gain. 

Options  to  improve  directional  control  may  be  categorized  into  expensive  and 
inexpensive  changes.  Inexpensive  options  include  those  changes  which  require  only  software  I- 
load  or  K-load  change  or  a  new  piloting  technique.  Software  I-loads  and  K-loads  are  constants 
which  are  not  coded  in  the  software  and  can  be  easily  changed  with  no  actual  software  code 
change.  Expensive  options  include  hardware  changes,  software  changes  which  require  a  code 
change,  and  changes  which  require  FCS  recertification. 

Table  5  summarizes  the  problems  with  the  current  FCS,  potential  changes  to  the  FCS  to 
reduce  susceptibility  to  PIO,  the  relative  cost  of  the  change,  and  whether  the  change  is  evaluated  in 
this  thesis. 

Table  5.  Potential  Changes  to  the  FCS 


Problem 

Potential  Fixes 

Relative  Cost 

Evaluated 

Slow  rudder  rate 

Hardware  modification 

Very  expensive 

No 

Slow  rudder  response 

Decrease  time  delays 

Linearize  rudder  command  shaper 

Disable  Ny  feedback 

Inexpensive 
Inexpensive 
Expensive 

Yes 
Yes 
Yes 

Powerful  rudder 

Decrease  pedal  gain 

Inexpensive 

Yes 

Position  of  rudder 

Redesign  Orbiter 

Very  expensive 

No 
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A.  LINEARIZE  THE  RUDDER  PEDAL  COMMAND  SHAPER 

The  rudder  command  shaper  was  described  in  Section  III.A.l.b.  Recall  that  the  shaping 
schedule         is         parabolic         and         conforms         to         the         following         function: 

DRMS  =  (0.042)DRMAN2  +  (0.\3\)DRMAN . 

The  shaper  can  be  linearized  by  changing  0.042  to  0  and  0. 13 1 .  to  1 .0.  The  above 
equation  becomes:  DRMS  =  DRMAN  .  This  change  provides  a  one-to-one  response  up  to  the 
limits  of  ±  22.5°  which  may  result  in  a  more  predictable  rudder  response. 

B.  REDUCE  THE  RUDDER  PEDAL  GAIN 

The  rudder  pedal  gain  was  described  in  Section  in. A. I.e.  Reducing  the  gain  will  reduce 
the  g-command  authority  on  the  rudder  and  may  reduce  over-controlling  of  the  Obiter.  The 
current  gain  value  is  0.08574  for  a  touchdown  speed  of  207  KEAS.  The  proposal  is  to  limit  the 
gain,  GNYDRM  to  0.04  at  683  fps  TAS  and  below.  This  change  will  reduce  the  gain  by  47% 
which  may  reduce  the  over-controlling  tendency. 

C.  DECREASE  TIME  LAGS 

The  time  lag  between  rudder  pedal  input  and  rudder  deflection  can  be  reduced  by 
changing  the  Ny  feedback  time  constant,  described  in  Section  III.  A.  1  f,  from  x  =  0.8  s  to  t  =  0.2  s. 
This  change  should  result  in  quicker  rudder  response  from  the  FCS  which  should  reduce  the 
susceptibility  to  PIO. 

D.  DISABLE  LATERAL  ACCELERATION  FEEDBACK 

Another  option  is  to  disable  Ny  feedback  in  combination  with  changing  the  Ny  feedback 
time  constant  from  t=0.8  s  to  t=0.2  s.  This  change  would  result  in  a  more  direct  rudder  FCS  by 
allowing  the  rudder  command  to  pass  directly  through  the  system  with  minimal  delay.  This 
change  would  require  software  code  change  and  is  relatively  expensive. 
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VI.  YAW  FCS  SIMULATION 

Directional  control  was  evaluated  by  incorporating  Orbiter  lateral-directional  control  laws 
and  longitudinal  and  lateral-directional  EOM  in  MATLAB  5.0  and  Simulink.  The  script  program, 
'eom.m',  defines  the  initial  conditions,  mass  properties,  stability  and  control  derivatives,  and 
longitudinal  and  lateral -directional  EOM  in  state  space  form.  The  code  for  'eom.m'  is  contained  in 
Appendix  A.  The  code  for  the  Simulink  simulation  is  contained  in  Appendix  B.  The  plotting 
program,  'FCSplot.m',  is  contained  in  Appendix  C. 

A.         SIMULATION  INITIAL  CONDITIONS 

The  case  investigated  is  comparable  to  the  STS-62  touchdown  flight  conditions.  The 
commander  input  approximately  +15°  (right)  rudder  pedal  deflection  followed  by  -15°  (left) 
rudder  pedal  deflection  in  the  two  point  stance  in  an  attempt  to  capture  runway  centerline. 
Landing  speed  was  approximately  207  KEAS  with  minimal  crosswind  (<  2  knots).  Table  6 
contains  the  initial  conditions  used  in  the  simulation. 


Table  6.  Initial  Conditions 

Parameter 

Initial  Value 

Vehicle  Weight 

228237  lbs 

Xcg 

1084.15  inches 

Ycg 

0 

Zcg 

370.65  inches 

Ixx 

8.58677e+05  slugs***2 

Iyy 

6.8699075e+06  slugs***2 

Izz 

7.1334485e+06slugs*ft2 

Wing  Surface  Area 

2690  ft2 

Chord  Length 

39.56  ft 

Ul 

207  KEAS  (353  fps) 

Runway  rolling  |x 

0.03 

Sea  Level  Density 

0.002377  slugs/ft3 

©i 

9° 

a 

10° 
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B. 


FCS  CONFIGURATIONS  EVALUATED 

Simulations  were  run  to  compare  six  FCS  configurations  as  defined  in  Table  7. 
Table  7.  FCS  Configurations  Evaluated 


Case  Number 

FCS  Configuration  (Simulation  Name) 

0 

Current  configuration  (baseline) 

1 

Disable  Ny  feedback  (nony) 

2 

Decrease  time  lag  t  to  0.2  (tau) 

3 

Linearize  rudder  shaper  (linrud) 

4 

Reduce  rudder  pedal  gain  (rudgain) 

5 

Cases  3,4,5  combined  (done) 

c. 


SIMULATION 


The  simulations  were  run  in  the  time  domain  for  30  seconds  using  Dormand-Prince 
integrator.  The  simulated  rudder  pedal  deflection  included  a  rudder  doublet  of + 1 5°  rudder  at  2 
seconds  followed  by  -15°  as  depicted  in  Figure  19.  The  master  Simulink  block  diagram  is  shown 
in  Figures  20. 
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Figure  19.  Rudder  Pedal  Deflection 
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Figure  20.  Simulink  Block  Diagram  of  Rudder  Doublet  Effect  on  Directional  Control 

The  lateral-directional  FCS  is  a  subsystem  of  the  entire  system  shown  in  Figure  20.  The 
FCS  is  labeled  'YAW  FCS*  in  Figure  2 1 .  This  subsystem  contains  the  simulation  code  and  block 
diagrams  for  the  FCS.  This  block  diagram  and  its  internal  block  diagrams  are  shown  in  Figure  2 1 . 


43 


-O-i 


Rudder 


< 

1 

« — 

353 
TAS 

GNY  DRM  Comp 

GNYDRM  COMP 


1 

0.8S+1 

3  RAY 

NYCOMP 

GERYFWDS 


ZxL     ~ *  GAIN1  Sum2       i> 


d0Ub,etRUDSHAPE 

A 


0    _l 
^pj    Ny  feedback 


0J 


Sum1 


f(u) 


►^> 


DRman(0.42|DRman|+0.131)       DRLr^r 
RUDSHAPE 


► 


RUDCMD 


NYCOMP 


Yaw  rate 


Yaw  Rate  Filter 


.7718479  -  0.3633895Z  _1  +0.49806552 
1-.3633895Z       +.26991 35z"2 


Yaw  Rate  Filter 


Figure  2 1 .  YAW  FCS  Block  Diagram 
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VII.  RESULTS 

Traces  of  Spedai,  8r,  Ny,  3,  r,  p,  \\i  and  <}>  versus  rime  are  shown  in  Figures  22-28.  All  six 
configurations  are  shown  on  each  plot  for  comparison  of  the  proposed  changes.  All  Simulink 
traces  shown  in  Figures  22-28  are  in  response  to  the  15°  doublet.  The  Simulink  model  was 
validated  by  comparing  the  results  with  actual  Shuttle  flight  data  and  Shuttle  Entry  Simulator 
(SES)  simulation  results. 

A.         MODEL  VALIDATION 

To  establish  validity  of  the  Simulink  model,  actual  STS-62  flight  data  and  SES  simulation 
results  were  compared  to  the  Simulink  results.  The  STS-62  flight  data  and  SES  results  are 
contained  in  Appendix  D.  The  STS-62  data  is  compared  to  the  Simulink  'Baseline'  traces,  since 
the  onboard  FCS  is  in  the  same  configuration.  The  STS-62  CDR  input  a  comparable  15°  rudder 
pedal  doublet  in  the  two  point  stance.  An  SES  evaluation,  using  the  same  initial  conditions  as  this 
simulation  and  the  STS-62  flight  data,  was  run  for  comparison.  Two  SES  runs  were  made.  The 
SES  'Old  FCS'  traces  contain  the  same  FCS  configuration  as  the  Simulink  'Baseline'  run.  The  SES 
'CF  FCS'  (Combined  Fix)  traces  contain  the  same  FCS  changes  as  the  Simulink  'Optimum  fix'. 

Rudder  deflection  traces  for  all  three  data  sets  (Simulink,  STS-62  flight  data,  and  SES) 
show  similar  magnitude  of  nearly  15°  deflection,  similar  response  time,  and  deflection  rate.  Ny 
traces  exhibit  non-minimal  phase  response  as  indicated  by  the  initial  low  magnitude  acceleration 
in  the  negative  direction.  Lag  times  are  similar  and  the  magnitude  of  the  Ny  traces  are 
comparable.  A  comparison  of  the  sideslip  angle  traces  indicate  the  same  initial  order-of- 
magnitude  3  resulting  from  the  doublet.  However,  the  Simulink  trace  shows  an  apparent  long 
term  directional  instability  which  may  be  the  reason  why  there  is  no  3  response  in  the  opposite 
direction.  The  long  term  directional  instability  is  reported  to  exist  in  the  actual  Orbiter  in  the  two 
point  stance.  This  instability  is  most  likely  not  as  strong  as  indicated  in  the  Simulink  model. 

Simulink  yaw  rate  and  roll  rate  traces  also  indicate  the  same  magnitude  and  lag  time  as 
SES  results  and  actual  flight  data.  The  SES  'Old  FCS'  roll  rate  data  is  'noisy'  due  to  the  gear 
model. 

This  comparison  of  flight  data  and  SES  results  with  these  Simulink  traces  results  in  a  high 
level  of  confidence  that  this  model  is  sufficient  to  analyze  and  compare  the  proposed  FCS 
configurations  for  short  term  Orbiter  directional  response. 
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B.         RUDDER  DEFLECTION 

The  rudder  deflection  plot  in  Figure  22  indicates  that  decreasing  the  time  constant  results 
in  the  quickest  rudder  response.  Simply  reducing  x  to  0,  as  in  the  'No  Ny  feedback'  option,  yields 
a  rudder  response  very  similar  to  the  pedal  command.  This  option  appears  the  best  with  respect  to 
rudder  surface  control,  yet  the  absence  of  Ny  feedback  generates  other  long  term  anomalies  in 
yaw  rate  and  siseslip  stability. 

Reducing  t  to  0.2  with  Ny  feedback,  as  in  the  'Shorter  time  constant'  trace,  yields  an 
acceptable  rudder  surface  response  with  the  benefits  of  Ny  feedback.  This  trace  shows 
approximately  one  half  second  of  time  delay  by  the  end  of  the  four  second  doublet.  This  time 
delay  could  also  be  considered  as  -45°  of  phase  lag.  While  45°  of  phase  lag  is  not  desired,  this 
option  appears  to  be  much  better  than  some  of  the  other  options.  The  baseline  configuration 
shows  nearly  two  seconds  of  time  delay  before  the  rudder  returns  through  the  zero  deflection 
position.  This  delay  equates  to  nearly  180°  of  phase  lag  by  the  end  of  the  maneuver.  This 
response  is  obviously  an  unacceptable  control  response  to  a  low  frequency  command. 

Linearizing  the  rudder  pedal  shaper,  with  no  additional  changes,  would  result  in  an 
increase  in  rudder  deflection  of  25%  over  the  baseline.  This  larger  rudder  deflection  for  the 
relatively  small  rudder  pedal  command  would  aggravate  the  already  difficult  fine  directional 
control.  Linearizing  the  rudder  with  no  other  changes  is  not  recommended. 

Reducing  the  rudder  pedal  gain,  with  no  additional  changes,  is  effective  in  decreasing  the 
magnitude  of  the  rudder  surface  response.  However,  the  rise  time  in  the  deflection  is  too  slow. 
The  slow  rise  time  distorts  the  predictability  of  the  rudder  response.  This  change  does  not  effect 
the  excessive  phase  lag. 

The  'optimum  fix',  which  includes  the  x=0.2  time  constant,  linearized  rudder  shaper.  and 
reduced  rudder  gain,  result  in  the  best  shaped  rudder  response  while  retaining  the  Ny  feedback. 
The  deflection  is  shaped  more  like  the  rudder  pedal  input  than  all  other  Ny  feedback  cases. 
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Rjdder  Deflection 
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Figure  22.  Rudder  Surface  Deflection 


C.         LATERAL  ACCELERATION  AT  THE  COMMANDERS  STATION 

The  lateral  acceleration  plot  shown  in  Figure  23  shows  that  eliminating  Ny  feedback 
results  in  greater  Ny  gain  and  more  time  required  to  return  to  a  zero  Ny  condition.  This  option 
displays  the  worst  long  term  directional  stability  characteristics. 

Both  the  'Shorter  time  constant'  and  'Linearized  rudder'  fixes  also  result  in  an  increase  in  g 
by  25%  over  the  baseline  FCS.   This  increase  in  lateral  acceleration  control  power  could  lead  to 
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directional  oversensitivity.  The  'Shorter  time  constant"  fix  results  in  the  least  phase  lag,  while 
'Linearized  rudder'  yields  the  longest  time  delay  of  the  Ny  feedback  options.  'Reduced  rudder 
gain'  results  in  an  appropriate  decrease  in  Ny  gain.  However,  this  change  alone  results  in 
excessive  time  delay. 

The  'Optimum  fix'  has  the  least  time  lag  and  acceptable  Ny  gain  characteristics. 
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Figure  23.  Lateral  Acceleration  at  the  Commander's  Station 


48 


D.         SIDESLIP 

The  sideslip  plot  shown  in  Figure  24  shows  the  total  rudder  control  power  in  the  classic 
sense  of  degrees  of  sideslip  per  degree  of  rudder  pedal.  Most  modern  aircraft  exhibit  rudder 
control  power  on  the  order  of  1 :4,  or  one  degree  of  sideslip  for  four  degrees  of  rudder  pedal.  The 
'Baseline'  Orbiter  trace  depicts  rudder  control  power  of  1:2.3.  This  apparent  oversensitivity  can 
lead  to  difficulty  in  precisely  controlling  sideslip. 

The  'No  Ny  feedback'  option  exhibits  excessive  rudder  control  power  with  an  additional 
long  term  instability.  This  long  term  response  is  unacceptable.  The  change  that  optimizes  the 
classic  1 :4  rudder  control  power  ratio  is  the  'Optimum  fix*. 

The  sideslip  plot  also  reveals  the  inherent  long  term  directional  instability  of  the  Orbiter  in 
the  two  point  stance.  This  instability  is  due  to  the  change  in  the  center  of  rotation  from  the  CG  to 
the  main  wheels  as  discussed  earlier. 
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Sdeslip 


Figure  24.  Sideslip  Angle 

E.        YAW  RATE  AND  HEADING 

Yaw  rate  and  heading  traces  are  shown  in  Figures  25  and  26.  Heading  was  derived  by 
integrating  yaw  rate.  Long  term  directional  instability  is  evident  in  both  yaw  rate  and  heading 
traces.  The  'Shorter  time  constant'  and  'Optimum  fix'  options  indicate  the  least  phase  lag  of  all 
options  with  Ny  feedback. 
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Figure  25.  Yaw  Rate 
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Figure  26.  Heading  Angle 


F. 


ROLL  RATE  AND  ANGLE  OF  BANK 


Adverse  roll  is  evident  as  shown  in  Figures  27  and  28.  This  response  is  undesirable 
because  adverse  roll  contributes  to  the  apparent  reduced  directional  control.  The  'Reduced  rudder 
gain'  results  in  the  least  adverse  roll.  However,  this  option  results  in  significant  phase  lag.  The 
"Optimum  fix'  offers  reduce  adverse  roll  and  the  least  phase  lag.   The  long  term  positive  dihedral 
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is  evident  in  the  long  term  positive  roll  rate.  This  roll  rate  is  in  response  to  the  long  term  negative 
sideslip. 
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Figure  27.  Roll  Rate 
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Figure  28.  Angle  of  Bank 
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G.         NATURAL  FREQUENCY  AND  DAMPING  RATIO 

The  so  called  'directional  natural  frequency',  tuD,  and  damping  ratio,  £d,  were  estimated 
from  the  plots  after  completion  of  the  doublet.  The  damping  ratio  was  estimated  using  the 
following  equation: 
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In 


<ZD  =  sin 


tan 


n 


Estimations  from  the  plots  indicate  that  mD=  0.6  rad/sec  and  £D  =  0.45.  As  expected,  the 
natural  frequency  and  damping  ratio  are  low.  This  low  natural  frequency  and  low  damping  ratio 
can  result  in  the  tendency  of  the  FCS  to  be  susceptible  to  PIO. 
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VIII.  CONCLUSIONS 

Orbiter  landings  indicate  both  long  term  directional  instability  and  the  potential  for  pilot 
induced  oscillations  during  landing  and  rollout  before  nosewheel  touchdown  with  the  baseline 
(current)  FCS.  The  Orbiter's  FCS  requires  improvements  to  increase  directional  control  in  the 
two  point  stance.  This  thesis  described  the  control  deficiency,  potential  improvements  to  the  FCS, 
and  evaluated  six  FCS  configurations  to  determine  their  effect  on  directional  control.  The 
evaluation  was  performed  by  modeling  the  Orbiter's  postlanding  lateral-directional  control  laws 
using  MATLAB  5.0  and  Simulink.  The  longitudinal  equations  of  motion  were  include  so  that  this 
model  can  more  readily  be  adapted  to  evaluate  longitudinal  control  problems. 

A.         RECOMMENDATIONS 

Of  the  six  FCS  configurations  evaluated,  the  option  which  provides  the  most 
improvement  over  the  baseline  configuration  is  the  'Optimum  fix'.  This  fix  includes  reducing  the 
time  lag  t  from  0.8  sec  to  0.2  sec,  linearizing  the  rudder  pedal  shaper,  and  reducing  the  rudder 
gain.  This  change  minimizes  the  apparent  phase  lag,  minimizes  the  adverse  roll,  reduces 
directional  oversensitivity,  and  provides  the  most  traditional  rudder  control  power.  This  option  is 
relatively  inexpensive  since  only  software  K-load  and  I-load  changes  are  required.  Software  code 
and  hardware  modifications  are  not  required.  The  data  indicates  that  PIO  tendencies  would  be 
reduced  over  the  baseline  system  if  this  option  were  implemented 

This  simulation  was  limited  to  analyzing  the  response  of  the  FCS  to  a  rudder  doublet 
since  the  initial  Orbiter  directional  difficulties  were  in  response  to  a  CDR  doublet.  The  doublet 
input,  in  essence,  is  a  frequency  domain  input.  Therefore,  frequency  domain  analysis  would  be 
appropriate.  However,  the  validity  of  the  Simulink  model  can  more  accurately  be  determined  by 
comparing  the  time  domain  response  traces  of  the  individual  states  to  actual  flight  data.  Although 
this  method  was  used  as  a  diagnostic  tool  to  study  a  frequency-concerned  phenomenon,  phase  lag 
can  be  estimated  from  the  traces.  In  addition,  these  changes  can  be  comapared  and  evaluated  with 
a  pilot  in  the  loop  in  time-based  simulators  such  as  the  SES  and  Ames  VMS. 

This  author  recommends  further  investigation  of  the  'Optimum  fix'  in  the  SES  and  Ames 
VMS.  Frequency  domain  analysis  would  also  be  appropriate.  In  addition,  the 
MATLAB/Simulink  program  used  for  this  evaluation  should  be  upgraded  to  include  a  gear  model, 
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drag  chute  model,  and  wind  model.     The  roll  and  pitch  channels  should  be  incorporated  for 
completeness. 

B.         USE  OF  COTS  SOFTWARE  FOR  ENGINEERING  EVALUATION 

This  author  concludes  that  commercially  available  engineering  software  packages  such  as 
MATLAB  may  be  used  to  analysis  proposed  FCS  changes  prior  to  full  up  simulation  in  more 
expensive  simulators.  Future  programs,  such  as  X38,  as  well  as  Orbiter  upgrade  projects  can 
benefit  from  the  use  of  PC-based  simulation. 
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APPENDIX  A.  MATLAB  CODE  FOR  SCRIPT  PROGRAM 

This  appendix  contains  the  Matlab  code  for  the  script  program,  "eom.m".  EOM.M 
defines  the  initial  conditions,  mass  properties,  stability  and  control  derivatives,  and  longitudinal 
and  lateral/directional  EOM  in  state  space  form. 

^Initial  conditions  and  EOM 

clear 

S  =  2690;    Isurface  area  of  the  wing  in  ft ^2 

C  =39.56;    %chord  length  in  ft 

B=78.057;    %wing  span  in  feet 

W=228237;    %weight  in  lbs 

g=32.174;    %gravitational  constant 

M=W/g;       %mass  in  slugs 

Ixx=858677; 

Iyy=6869907.5; 

Izz=7133448.5; 

%Shift  the  moments  of  inertia  from  the  eg  to  center  of  rotation 
% (between  the  main  gear) 
Izmg=Izz+M* ((1172. 6-1084. 15) /12) "2 

%b/h=0.05  at  alpha  =  10 

Mu=0.03;  %runway  rolling  Mu 

Ul=353;  %initial  velocity  in  the  x  direction  ft/sec 

rho=  0.002377;  %sea  level  pressure 

theta=9*pi/180  %Pitch  angle  in  radians 

alpha=10*pi/180  %Angle  of  attack  in  radians 

q=0 .5*rho*Ul~2  %Dynamic  pressure  at  sea  level 

%Non-Dimensional  Stability  and  Control  Derivatives 

CLa=0. 05948; 

CDa=0. 009924 ; 

CL1=0. 54055; 

CD1=0. 16712; 

CLu=0; 

CDu=0; 

Cyb=-0. 01705; 
Cyp=0; 
Cyr=0; 

Cydeltar=0. 00138; 
Clb=-0. 00205; 
Clp=-0.3189; 
Clr=0.1890; 
Cldeltar=0. 00032; 
Cnb=0.0010; 
Cnp=0.028; 
Cnr=-0.2  925; 
Cndeltar=- 0.00054; 

%Longitudinal  Dimensional  Stability  and  Control  Derivatives 
Xu=-q*S* (CDu+2*CDl) / (M*U1) 
Xa=-q*S* (CDa+2*CLl) /M 
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Xq=0; 

Xdeltae=0; 

Zu=-q*S* (CLu+2*CLl) / (M*U1 

Za=-q*S* (CLa+2*CDl) /M 

Zq=0; 

Zdeltae=0; 

Mu=0; 

Ma=0, 

Mq=0, 

Mdeltae=0; 

%Lateral-directional  Dimensional  Stability  and  Control  Derivatives 

Yb=q*S*Cyb/M 

Yp=q*S*B*Cyp/ (2*M*U1)    Sis  zero 

Yr=q*S*B*Cyr/ (2*M*U1)    Sis  zero 

Ydeltaa=0; 

Ydeltar=q*S*Cydeltar/M 

Lb=q*S*B*Clb/Ixx 
Lp=q*S*B~2*Clp/ (2*Ixx*Ul) 
Lr=q*S*BA2*Clr/ (2*Ixx*Ul) 
Ldeltaa=0; 
Ldeltar=q*S*B*Cldeltar/Ixx 

Nb=q*S*B+Cnb/Izmg 
Np=q*S*B"2*Cnp/ (2*Izmg*Ul) 
Nr=q*S*B"2*Cnr/ (2*Izmg*Ul) 
Ndeltaa=0; 
Ndeltar=q*S*B*Cndeltar/Izmg 

SDefinition  of  dx/dt=  Ax+Bu 

%xdot=[udot  alphadot  qdot  thetadot  betadot  pdot  rdot  phidot] 
Sx=[u  alpha  q  theta  beta  p  r  phi] 
%U=[DELTAE  DELTAA  DELTAR] 

A  =  [Xu  Xa  Xq  -g*cos (theta)  0  0  0  0; 

Zu/Ul  Za/Ul  Zq/Ul+1  -g*sin (theta) /Ul  0  0  0  0; 

Mu  Ma  Mq  0  0  0  0  0; 

00100000; 

0  0  0  0  Yb/Ul  Yp/Ul  Yr/Ul-1  g*cos (theta) /Ul ; 

0  0  0  0  Lb  Lp  Lr  0; 

0  0  0  0  Nb  Np  Nr  0; 

0  0  0  0  0  1  tan (theta)  0] 

B=[Xdeltae  0  0; 

Zdeltae/Ul  0  0; 

Mdeltae  0  0; 

0  0  0; 

0  Ydeltaa/Ul  Ydeltar/Ul; 

0  Ldeltaa  Ldeltar; 

0  Ndeltaa  Ndeltar; 

0  0  0] 

%  dy/dt  =  Cx+Du;  C=  identity  matrix  D=0  so  that  the  output 
%of  the  state  space  is  X 
C=[eye(8)  ] 
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D=[zeros (8,3) ] 

B=[Xdeltae  0  0; 

Zdeltae/Ul  0  0; 

Mdeltae  0  0; 

0  0  0; 

0  Ydeltaa/Ul  Ydeltar/Ul; 

0  Ldeltaa  Ldeltar; 

0  Ndeltaa  Ndeltar; 

0  0  0] 
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APPENDIX  B.  MATLAB  CODE  FOR  SIMULATION 

The  code  for  the  actual  simulation,  "baseline.m"  is  contained  in  this  Appendix.  The  code 
is  for  the  baseline  system. 


Model  { 

Name 

"baseline" 

Version 

200 

SimParamPage 

SampIeTimeCoIors 

WideVectorLines 

Solver 

off 

off 

PaperOrientation 
StartTime 

landscape 
"0.0" 

StopTime 

Solver                  ode45 

"30" 

RelTol                  "01" 

AbsTol                 "le-6" 

Refine                  "1" 

MaxStep 
InitialStep 
Fixed  Step 
MaxOrder 

"05" 

"auto" 

"auto" 

5 

OutputOption 
OutputTimes 
LoadExtemallnput 
External  Input 
SaveTime 

RefineOutputTimes 

off 

"ft  u]" 
off 

TimeSaveName 

"tout" 

SaveState 

off 

StateSaveName 

"xout" 

SaveOutput 
OutputSaveName 
LoadlnitialState    off 

off 

"yout" 

InitialState 

"xlnitial" 

SaveFinal  State      off 

FinalStateName    "xFinal' 

t 

LimitMaxRows 

off 

MaxRows 

"1000" 

Decimation 

1*  1  11 

AlgebraicLoopMsg 

MinStepSizeMsg 

UnconnectedlnputMsg 

UnconnectedOutputMsg 

UnconnectedLineMsg 

ConsistencyChecking 

ZeroCross 

warning 

warning 

warning 

warning 

warning 

off 

on 

BlockDefaults  { 
Orientation 
ForegroundColor 

right 
black 
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BackgroundColor 

white 

DropShadow 

off 

NamePlacement 

normal 

FontName 

"Helvetica" 

FontSize 

10 

FontWeight 

normal 

FontAngle 

normal 

ShowName 

> 

on 

AnnotationDefaults  { 

Horizontal  Alignment 

center 

VerticalAlignment 

middle 

ForegroundColor 

black 

BackgroundColor 

white 

DropShadow 

off 

FontName 

"Helvetica" 

FontSize 

10 

FontWeight 

normal 

FontAngle 

normal 

System  { 

Name                  "baseline" 

Location 

[0,  3,  644,  447] 

Open                   on 

ScreenColor 

white 

Block  { 

BlockType 

Derivative 

Name 

"BETADOT" 

Position 

[385,320,415,340] 

Orientation           left 

} 
Block  { 

BlockType 

Clock 

Name 

"Clock" 

Position 

[30,  330,  50,  350] 

DeleteFcn 

"simclock  BlocklsBeingDestroyed" 

PostSaveFcn         "simclock  Save" 

Location 

[12,31,  119,76] 

i 
Block  { 

BlockType 

Demux 

Name 

"Demux" 

Position 

[370,  154,405,236] 

Outputs 

"8" 

i 
Block  { 

BlockType 

Mux 

Name 

"Mux" 

Position 

[310,253,335,347] 

Orientation           left 

Inputs 

"3" 
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Block  { 

BlockType 

Mux 

Name 

"Muxl" 

Position 

[235,  178,270,212] 

Inputs 

„3„ 

i 
Block  { 

BlockType 

Scope 

Name 

"Ny" 

Position 

[215,395,245,425] 

Orientation 

down 

NamePlacement 

alternate 

Floating 

off 

Location 

[188,235,512,474] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

»-5» 

YMax 

"5" 

SaveTo  Workspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 
} 

"0" 

Block  { 
BlockType 

Fen 

Name 

"Ny  " 

Position 

[240,  277,  275,  323] 

Orientation 

left 

Expr 

"-(353 

*u[3]+353*u[2]+u[l]*< 

Block  { 

BlockType 

Integrator 

Name 

"PSI" 

Position 

[505,  380,  525,  400] 

ExternalReset 

none 

InitialConditionSource 

internal 

InitialCondition 

"0" 

LimitOutput 

off 

UpperSaturationLimit 

"inf" 

LowerSaturationLimit 

"-inf 

ShowSaturationPort 

off 

ShowStatePort 

off 

AbsoluteTolerance 

} 

"auto" 
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Block  { 
BlockType 
Name 
Position 

ShowPortLabels 
MaskType 
MaskDescription 
MaskHelp 


MaskPrompt  String 

MaskStyleString 

MaskVariables 

Masklnitialization 


"\n" 


Subsystem 

"PuIseVnGeneratorl " 

[25,  35,  55,  65] 

on 

"Pulse  Generator" 

"Pulse  Generator" 

"Pulse  generator  which  ensures  pulse  transitions  " 

"are  hit.  Provides  a  vector  of  pulses  when  " 

"the  height  is  entered  as  a  vector.  Unmask  " 

"to  see  how  it  works. " 

"Period  (secs):|Duty  cycle  (%  of  period):  |Amplitu" 

"de:|Start  time:" 

"edit,edit,edit,edit" 

"period=@  1  ;duty=@2,arnplitude=@3,start=@4," 


MaskDisplay         "plot(0,0,100,100,[95,90,90,75,75,60,60,45,45,30," 

"30, 1 5, 1 5,5],[25,25,75,75,25,25,75,75,25,25,75,75' 
",25,25]);  \n 


\n 


MasklconFrame 

on 

MasklconOpaque 

on 

MasklconRotate 

off 

MasklconUnits 

autoscale 

MaskValueString 

"1000|.2|15|5" 

System  { 

Name 

"Pulse\nGeneratorl" 

Location 

[205,  368,  544,  559] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Constant 

Name 

"Constant  1" 

Position 

[40,  135,60,  155] 

Value 

> 

"0" 

J 
Block  { 

BlockType 

Sin 

Name 

"Sine  Wave" 

Position 

[35,  83,  65,  107] 

Amplitude 

ff  i  ii 

Frequency 

"2*pi./period" 

Phase 

"-2*pi*(start  /period  +  duty/200  -  1/4)" 

SampleTime 

»_,„ 

Block  { 

BlockType 

Step 

Name 

"Step" 
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Position 

Time 

Before  "0" 

After 

} 
Block  { 

BlockType 

Name 

Position 

Threshold 

} 
Block  { 

BlockType 

Name 

Position 

Port 

Output  WhenDisabled 

InitialOutput 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 
Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 
Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 
Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 
} 


[35,  27,  65,  53] 
"start" 

"amplitude" 


Switch 
"Switch" 

[120,65,  155,  125] 
"sin(pi*(0.5  -  duty/100))" 


Out  port 

"outl" 

[250,88,280,  102] 

"1" 

held 

"0" 


} 
} 
Block  { 

BlockType 

Name 


"Sine  Wave" 
1 
[5,  0;  20,  0] 

"Switch" 

2 


"Switch" 


"outl" 
1 


"Constant  1" 


[25,  0;  0,  -30] 

"Switch" 
3 


"Step" 


1 
[20,  0;  0,  35] 

"Switch" 
1 


Subsystem 
"Pulse\nGenerator2" 
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Position 

ShowPortLabels 

MaskType 

MaskDescription 

MaskHelp 


MaskPromptString 


[25,  95,  55,  125] 

on 

"Pulse  Generator" 

"Pulse  Generator" 

"Pulse  generator  which  ensures  pulse  transitions  " 

"are  hit.  Provides  a  vector  of  pulses  when  " 

"the  height  is  entered  as  a  vector.  Unmask  " 

"to  see  how  it  works." 

"Period  (secs):|Duty  cycle  (%  of  period):  |Amplitu" 

"de: | Start  time:" 

"edit,edit,edit,edit" 

"period=@l;duty=@2;amphtude=@3;start=@4;" 


MaskStyleString 

MaskVariables 

Masklnitialization       "\n" 

MaskDisplay         "plot(0,0,100,100,[95,90,90,75,75,60,60,45,45,30," 

"30, 1 5, 1 5,5],[25,25,75,75,25,25,75,75,25,25,75,75" 

",25,25]);  \n 


\n 


MasklconFrame 

on 

MasklconOpaque 

on 

MasklconRotate 

off 

MasklconUnits 

autoscale 

MaskValueString 

"1000|.2|-15|7" 

System  { 

Name 

"PuIse\nGenerator2" 

Location 

[205,  368,  544,  559] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Constant 

Name 

"Constant  1" 

Position 

[40,  135,60,  155] 

Value 

"0" 

} 
Block  { 

BlockType 

Sin 

Name 

"Sine  Wave" 

Position 

[35,  83,  65,  107] 

Amplitude 

"1" 

Frequency 

"2*pi./period" 

Phase 

"-2*pi*(start./period  +  duty/200  -  1/4)" 

SampleTime 

H       1  tl 

i 
Block  { 

BlockType 

Step 

Name 

"Step" 

Position 

[35,  27,  65,  53] 

Time 

"start" 

Before 

"0" 
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After 


"amplitude" 


Block  { 
BlockType 

Switch 

Name 

"Switch" 

Position 

[120,65,  155,  125] 

Threshold 

"sin(pi*(0.5  -  duty/1 00)) 

) 
Block  { 

BlockType 

Outport 

Name 

"outl" 

Position 

[250,88,280,  102] 

Port 

"1" 

OutputWhenDisabled 

held 

InitiaJOutput 

} 
Line  { 

"0" 

SrcBlock 

"Step" 

SrcPort 

1 

Points 

[20,  0: 

,  0,  35] 

DstBlock 

"Switch" 

DstPort 

1 

) 
Line  { 

SrcBlock 

"Constant  1" 

SrcPort 

1 

Points 

[25,  0 

,  0,  -30] 

DstBlock 

"Switch" 

DstPort 

} 
Line  { 

3 

SrcBlock 

"Switch" 

SrcPort 

1 

DstBlock 

"outl" 

DstPort 

} 
Line  { 

1 

SrcBlock 

"Sine  Wave" 

SrcPort 

1 

Points 

[5,0; 

20,0] 

DstBlock 

"Switch" 

DstPort 

2 

Block  { 
BlockType 

Derivative 

Name 

"RDOT" 

Position 

[385,260,415,280] 

Orientation 

} 

left 
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Block  { 

BlockType 

Scope 

Name 

"Scope" 

Position 

[395,  435,  425, 

465] 

Floating 

off 

Location 

[115,431,439, 

670] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-5" 

YMax 

iicn 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

tf  -i  it 

Samplelnput 

off 

SampleTime 

"0" 

Block  { 

BlockType 

Scope 

Name 

"Scopel" 

Position 

[390,  380,  420, 

410] 

Floating 

off 

Location 

[100,  161,424, 

400] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-5" 

YMax 

"5" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 

"0" 

i 
Block  { 

BlockType 

StateSpace 

Name 

"State-space" 

Position 

[295,  180,  355, 

210] 

A 

"A" 

B 

"B" 

C 

llfl|l 

D 

"D" 

XO 

"[0;l 

3;0;0;0;0,0;0]" 
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Block  { 

BlockType 

Sum 

Name 

"Sum" 

Position 

[100,77,  130,  108] 

Inputs 

"++" 

t 
Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace" 

Position 

[95,  145,  155,  175] 

VariableName 

"doublet" 

Buffer 

"inf" 

Decimation           "  1 

It 

SampleTime          "-3 

} 
Block  { 

:" 

BlockType 

ToWorkspace 

Name 

"To  Workspace  1" 

Position 

[175,250,235,280] 

VariableName 

"dr" 

Buffer 

"inf" 

Decimation           "  1 

n 

SampleTime         "-] 
\ 

[" 

Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace2" 

Position 

[95,365,  155,395] 

Orientation           left 

VariableName 

"Ny" 

Buffer 

"inf' 

Decimation           "  1 

It 

SampleTime         "-] 

1" 

Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace3" 

Position 

[540,  110,600,  140] 

VariableName 

"beta" 

Buffer 

"inf 

Decimation           "  1 

n 

SampleTime          "-1 

I" 

i 
Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace4" 

Position 

[570,  150,630,  180] 

VariableName 

"p" 

Buffer 

"inf' 

Decimation           "  1 

It 
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SampleTime          "-1" 

X 

) 
Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace5 " 

Position 

[585,  220,  645,  250] 

VariableName 

Hj.ii 

Buffer 

"inf 

Decimation           "  1 " 

SampleTime         "-1" 

x 

Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace6" 

Position 

[585,  335,  645,  365] 

VariableName 

"psi" 

Buffer 

"inf 

Decimation           "  1 " 

SampleTime          "- 1 " 

x 

Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace7" 

Position 

[70,  320,  130,  350] 

VariableName 

"time" 

Buffer 

"inf 

Decimation            "  1 " 

SampleTime         "-1" 

x 

Block  { 

BlockType 

ToWorkspace 

Name 

"To  Workspace8" 

Position 

[585,280,645,310] 

VariableName 

"aob" 

Buffer 

"inf 

Decimation           "1" 

SampleTime         "-1" 

X 

f 
Block  { 

BlockType 

Subsystem 

Name 

"YAW  FCS" 

Position 

[175,  179,205,231] 

ShowPortLabels 

off 

System  { 

Name 

"YAW  FCS" 

Location 

[4,  44,  648,  440] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Inport 

Name 

"in  1" 
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Position 

[35,  180,  50,  200] 

Port 

ll  1  II 

Port  Width 

„_j„ 

SampleTime 

.,_,., 

* 
Block  { 

BlockType 

Inport 

Name 

"in  2" 

Position 

[100,290,  120,310] 

Port 

"2" 

Port  Width 

ll     i  ll 

SampleTime 

»_1« 

i 
Block  { 

BlockType 

Inport 

Name 

"in_3" 

Position 

[290,270,310,290] 

Port 

H ->  il 

PortWidth 

ii     i  « 

SampleTime 

ll      |  ll 

Block  { 

BlockType 

Subsystem 

Name 

"COSALF" 

Position 

[320,  70,  350,  120] 

Orientation 

down 

ShowPortLabels 

off 

System  { 

Name 

"COSALF" 

Location 

[318,77,607,247] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

MATLABFcn 

Name 

"Cos(  Alpha)" 

Position 

[57,95,  133,  125] 

Orientation 

down 

MATLABFcn 

"cos" 

OutputWidth 
\ 

M„ 

i 
Block  { 

BlockType 

Constant 

Name 

"alpha" 

Position 

[58,50,  132,65] 

Orientation 

down 

Value 

"10*pi/180" 

} 
Block  { 

BlockType 

Outport 

Name 

"out  1" 

Position 

[85,  150,  105,  170] 
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Orientation  down 

Port  "1" 

Output  WhenDisabled 
InitialOutput         "0" 


held 


} 


Line  { 

SrcBIock 

"Cos(  Alpha)" 

SrcPort 

1 

DstBlock 

"out_l" 

DstPort 

} 
Line  { 

1 

SrcBIock 

"alpha" 

SrcPort 

1 

DstBlock 

"Cos(  Alpha)" 

DstPort 
} 
> 
} 
Block  { 

1 

BlockType 

Constant 

Name 

"Constant" 

Position 

[60,  125,  80,  145] 

Value 

> 

"04" 

Block  { 

BlockType 

Subsystem 

Name 

"DRERR" 

Position 

[335,  255,  365,  305] 

ShowPortLabels 

off 

System  { 

Name 

"DPvERR" 

Location 

[180,294,476,476] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Inport 

Name 

"in  1" 

Position 

[15,75,35,95] 

Port 

Ml  II 

PortWidth 

tf      I  M 

SampleTime 

II       1  II 

i 
Block  { 

BlockType 

DiscreteFilter 

Name 

"YAW  RATE  BBF" 

Position 

[60,56,245,  114] 

Numerator 

"[.7718479  -.3633895.4980655] 

Denominator 

"[1- 

.3633895  .2699135]" 

SampleTime 

} 

•i  1  ii 
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Block  { 
BlockType 
Name 
Position 

Port  "1" 

OutputWhenDisabled 
InitialOutput         "0" 


Outport 
"out_l" 
[265,  75,  285,  95] 

held 


Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 
} 
} 
} 

Block  { 
BlockType 
Name 
Position 
Value 

} 
Block  { 

BlockType 

Name 

Position 

Inputs 

} 
Block  { 

BlockType 

Name 

Position 

Inputs 

} 
Block  { 

BlockType 

Name 

Position 

Orientation 

Value 

} 
Block  { 

BlockType 

Name 

Position 


"YAWRATEBBF" 

1 
"out_l" 

1 


"in_l" 

1 

"YAW  RATE  BBF 

1 


Constant 

"DRT" 

[135,245,  150,265] 

"0" 


Product 
"GAIN" 
[135,  175,  160,  195] 


Product 

"GAIN" 

[355,  185,380,205] 


Constant 
"GDRC" 

[475,300,515,330] 

up 

"6.325" 


Constant 
"GERYFWDS" 

[485,  105,  505,  125] 
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Orientation 

down 

Value 
} 

'•.5" 

Block  { 
BlockType 

Subsystem 

Name 

"GNYDRM  COMP" 

Position 

[105,65,  135,  115] 

Orientation 

down 

ShowPortLabels 

off 

System  { 

Name 

"GNYD 

Location 

[355,  132,  620,  292] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Lookup 

Name 

"GNY  DRM  Comp" 

Position 

[75,43,  150,  107] 

Orientation 

left 

InputValues 

"[0 

250    900    10000]" 

Output  Values 
i 

"[0.1 

0.1   0.01   001]" 

Block  { 

BlockType 

Constant 

Name 

"TAS" 

Position 

[170,75,  190,95] 

Orientation 

left 

Value 

» 

"353" 

1 
Block  { 

BlockType 

Outport 

Name 

"GNYDRM" 

Position 

[15,70,35,90] 

Orientation 

left 

Port 

"1" 

OutputWhenDisabled 

held 

InitialOutput 

"0" 

t 
Line  { 

SrcBlock 

"GNY  DRM  Comp" 

SrcPort 

1 

DstBlock 

"GNYDRM" 

DstPort 

1 

Line  { 

SrcBlock 

"TAS" 

SrcPort 

1 

DstBlock 

"GNY  DRM  Comp" 

DstPort 

1 
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BlockType 

Subsystem 

Name 

"NYCOMP" 

Position 

[295,  170,325,220] 

ShowPortLabels 

off 

System  { 

Name 

"NYC01 

Location 

[227,  109,  594,  284] 

Open 

off 

ScreenCoIor 

white 

Block  { 

BlockType 

Inport 

Name 

"in  1" 

Position 

[15,80,35,  100] 

Port 

"1" 

PortWidth 

»_1" 

SampleTime 

It         1     H 

Block  { 

BlockType 

Gain 

Name 

"GRAY" 

Position 

[215,77,240,  103] 

Gain 

IIC« 

Block  { 

BlockType 

TransferFcn 

Name 

"Transfer  Fen" 

Position 

[80,65,  175,  115] 

Numerator 

"[1]" 

Denominator 

} 
Block  { 

"[0.8  1]" 

BlockType 

Outport 

Name 

"out  1" 

Position 

[260,  80,  280,  100] 

Port 

it  i  II 

Output  WhenDisabled 

held 

Initial  Output 

\ 

"0" 

f 
Line  { 

SrcBlock 

"in  1" 

SrcPort 

1 

DstBlock 

"Transfer  Fen" 

DstPort 

} 
Line  { 

1 

SrcBlock 

"GRAY" 

SrcPort 

1 

DstBlock 

"out  1" 

DstPort 

1 
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ine  { 
SrcBlock 

"Transfer  Fen" 

SrcPort 

1 

DstBlock 

"GRAY" 

DstPort 

1 

Block  { 
BlockType 

Product 

Name 

"Product  1" 

Position 

[510,  177,535,213] 

Inputs 

"3" 

Block  { 

BlockType 

Saturate 

Name 

"RUDCMD" 

Position 

[555,  190,580,210] 

UpperLimit 

"27.1" 

LowerLimit 

} 

"-27. 1" 

Block  { 
BlockType 

Subsystem 

Name 

"RUDSHAPE" 

Position 

[80,  165,  110,215] 

ShowPortLabels 

off 

System  { 

Name 

"RUDSHAPE" 

Location 

[75,  132,  474,  272] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Inport 

Name 

"in  J" 

Position 

[15,  70,  35,  90] 

Port 

"1" 

PortWidth 

„_,„ 

SampleTime 

} 
Block  { 

(l      1  " 

BlockType 

Saturate 

Name 

"DR  Limiter" 

Position 

[210,66,235,94] 

Font  Size 

8 

UpperLimit 

"+22.5" 

LowerLimit 

"-22.5" 

Block  { 

BlockType 

Fen 

Name 

"DRman(0  42|DRi 
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Position 

[105,70,  145,90] 

Expr 

"(abs(u)*0.042+.131)*u" 

J 
Block  { 

BlockType 

Outport 

Name 

"out  1" 

Position 

[265,  70,  285,  90] 

Port 

ii  |  H 

OutputWhenDisabled 

held 

InitialOutput 

"0" 

Line  | 

SrcBlock 

"DR  Limiter" 

SrcPort 

1 

DstBlock 

"out   1" 

DstPort 

> 

1 

Line  { 

SrcBlock 

"in J" 

SrcPort 

1 

DstBlock 

"DRman(0.42|DRman|+0. 131)" 

DstPort 

1 

Line  { 

SrcBlock 

"DRman(0.42|DRman|+0. 131)" 

SrcPort 

1 

DstBlock 

"DR  Limiter" 

DstPort 

} 
} 

1 

} 

Block  { 
BlockType 

Scope 

Name 

"Scope" 

Position 

[400,  130,430,  160] 

Floating 

off 

Location 

[188,365,512,604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

yonly 

TimeRange 

"auto" 

YMin 

"-5" 

YMax 

llfli 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

H  1  tl 

Samplelnput 

off 

SampleTime 

"0" 
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Block  { 

BlockType 

Scope 

Name 

"Scopel" 

Position 

[275,  110,305, 

140] 

Floating 

off 

Location 

[188,365,512, 

604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-5" 

YMax 

"5" 

SaveTo  Workspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

ft  i  ff 

Samplelnput 

off 

SampleTime 

} 
Block  { 

"0" 

BlockType 

Scope 

Name 

"Scope2" 

Position 

[425,  300,  455, 

330] 

Floating 

off 

Location 

[188,365,512, 

604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-5" 

YMax 

»r  tt 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

11  1  11 

Samplelnput 

off 

SampleTime 
1 

"0" 

Block  { 

BlockType 

Scope 

Name 

"Scope3" 

Position 

[215,  115,245, 

145] 

Floating 

off 

Location 

[188,365,  512, 

604] 

Open 

off 

Grid 

on 
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TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-5" 

YMax 

"5" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 

} 

"0" 

Block  { 
BlockType 

Scope 

Name 

"Scope4" 

Position 

[240,  250,  270, 

280] 

Floating 

off 

Location 

[188,365,512, 

604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-5" 

YMax 

"5" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 

"0" 
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Block  { 

BlockType 

Sum 

Name 

"Sum" 

Position 

[420,  190,440, 

210] 

Inputs                    "++" 

Block  { 

BlockType 

Sum 

Name 

"Suml" 

Position 

[235,  177,  255, 

208] 

Inputs                    "-+" 

Block  f 

BlockType 

Sum 

Name 

"Sum2" 

Position 

[180,  175,200, 

195] 

Inputs                    "++" 
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} 

Block  { 

BlockType 

Subsystem 

Name 

"YBB  FILTER" 

Position 

[445,  15,475,65] 

ShowPortLabels 

off 

System  { 

Name 

"YBB 

_FILTER" 

Location 

[55,  69,  652,  308] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Inport 

Name 

"in  1" 

Position 

[15,90,35,  110] 

Port 

ft -j  ft 

PortWidth 

"_1» 

SampleTime 

f      i  fl 

Block  { 

BlockType 

DiscreteFilter 

Name 

"FUterl" 

Position 

[60,76,265,  124] 

Numerator 

"[.699558  -.82766  +.5599439] 

Denominator 

"[1   - 

•1.08273  +.554071]" 

SampleTime 

"  i  It 

i 
Block  { 

BlockType 

DiscreteFilter 

Name 

"Filter2" 

Position 

[325,  79,  525,  121] 

Numerator 

"[.940726  -1.04423   .801291]" 

Denominator 

"[1   - 

96303  .66162]" 

SampleTime 

It  i  it 

i 
Block  { 

BlockType 

Outport 

Name 

"out   1" 

Position 

[550,90,  570,  110] 

Port 

ft  -I  ft 

OutputWhenDisabled 

held 

InitialOutput 

1 

"0" 

Line  { 

SrcBlock 

"Filter2" 

SrcPort 

1 

DstBlock 

"out  1" 

DstPort 

t 

1 

Line  { 

SrcBlock 

"in  1" 

82 


SrcPort 

1 

DstBlock 

"Filter  1" 

DstPort 

1 

Line  { 

SrcBlock 

"Filterl" 

SrcPort 

1 

DstBlock 

"Filter2" 

DstPort 

1 

} 

} 

} 

Block  { 

BlockType 

Outport 

Name 

"out   1" 

Position 

[605,  190,625,210] 

Port 

M  1   fl 

Output  WhenDisabled 

held 

InitialOutput 

} 
Line  { 

"0" 

SrcBlock 

"Sum" 

SrcPort 

1 

Points 

[50,  0] 

DstBlock 

"Product  1" 

DstPort 

} 

Line  { 

2 

SrcBlock 

"in  2" 

SrcPort 

1 

Points 

[95,  0] 

DstBlock 

"Suml" 

DstPort 

} 
Line  { 

2 

SrcBlock 

"in  1" 

SrcPort 

1 

DstBlock 

"RUDSHAPE" 

DstPort 

} 
Line  { 

1 

SrcBlock 

"DRT" 

SrcPort 

1 

Points 

[10,  0] 

DstBlock 

"Sum2" 

DstPort 

} 
Line  { 

2 

SrcBlock 

"RUDSHAPE" 

SrcPort 

1 

DstBlock 

"GAIN" 
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DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"GAIN" 
1 

"Sum2" 
1 


Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


"NYCOMP" 
1 

[5,  0;  0,  5] 

"GAIN" 

2 


Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


[0,  65] 


"COSALF" 
1 

"GAIN" 

1 


"Product  1" 

1 

"RUDCMD" 

1 


Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"GERYFWDS" 

1 
"Product  1" 

1 


"GDRC" 

1 

"Product  1" 
3 


"RUDCMD" 

1 

"out  1" 


Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 


[0,0] 


»n  a  txi    » 


'GAIN 
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DstBlock 

"Sum" 

DstPort 

< 

1 

1 
Branch  { 

DstBlock 

"Scope" 

DstPort 

1 

} 

} 

Line  { 

SrcBlock 

"Suml" 

SrcPort 

1 

Points 

[0,0] 

Branch  } 

DstBlock 

"NYCOMP" 

DstPort 

1 

Branch  { 

DstBlock 

"Scopel" 

DstPort 

1 

} 

} 

Line  { 

SrcBlock 

"in_3" 

SrcPort 

1 

DstBlock 

"DRERR" 

DstPort 

1 

} 

Line  { 

SrcBlock 

"DRERR" 

SrcPort 

1 

Points 

[20,  0] 

Branch  { 

Points 

[15,0] 

DstBlock 

"Sum" 

DstPort 

> 

2 
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Branch  { 

Points 

[0,  35] 

DstBlock 

"Scope2" 

DstPort 

} 

1 
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Line  { 

SrcBlock 

"Sum2" 

SrcPort 

1 

Points 

[0,0] 

Branch  { 

DstBlock 

"Suml" 

DstPort 

1 

} 

Branch  { 
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Points 
DstBlock 

[0,  -55] 

"Scope 

DstPort 

1 

} 
} 
Line  { 

SrcBlock 

"GNYD 

SrcPort 

1 

DstBlock 

"GAIN" 

DstPort 

1 

Block  { 

BlockType 

Scope 

Name 

"alfa" 

Position 

[415,25,445,55] 

Orientation 

up 

NamePIacement 

alternate 

Floating 

off 

Location 

[188,365,512,604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

yonly 

TimeRange 

"auto" 

YMin 

"0  15" 

YMax 

"0.25" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

M  1  II 

Samplelnput 

off 

SampleTime 

X 

"0" 

s 
Block  { 

BlockType 

Scope 

Name 

"aob" 

Position 

[550,  300,  580,  330] 

Floating 

off 

Location 

[188,235,512,474] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-01" 

YMax 

"0.02" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 
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LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

tl   1   M 

Samplelnput 

off 

SampleTime 
} 

"0" 

Block  { 
BlockType 

Scope 

Name 

"beta" 

Position 

[505,  155,535,  185] 

Floating 

off 

Location 

[188,235,512,474] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-0.4" 

YMax 

"0" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 

} 
Block  { 

"0" 

BlockType 

Constant 

Name 

"deltaa" 

Position 

[190,  135,210,  155] 

Value 

"0" 

i 
Block  { 

BlockType 

Constant 

Name 

"deltae" 

Position 

[190,95,210,  115] 

Value 

} 
Block  { 

"0" 

BlockType 

Scope 

Name 

"doublet" 

Position 

[145,30,  175,60] 

Orientation 

up 

NamePlacement 

alternate 

Floating 

off 

Location 

[188,365,512,604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 
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TimeRange 

"auto" 

YMin 

"-5" 

YMax 

"5" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

ii  i  it 

Samplelnput 

off 

SampleTime 

"0" 

i 
Block  { 

BlockType 

Scope 

Name 

"dr" 

Position 

[240,  230,  270, 

260] 

Floating 

off 

Location 

[188,365,  512, 

604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

yonly 

TimeRange 

"auto" 

YMin 

It     CM 

YMax 

"5" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 

\ 

"0" 

) 
Block  { 

BlockType 

Scope 

Name 

"psi" 

Position 

[565,  375,  595, 

405] 

Floating 

off 

Location 

[188,235,512, 

474] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-0.05" 

YMax 

"0.15" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 
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SampleTime 

X 

"0" 

i 
Block  { 

BlockType 

Scope 

Name 

"q" 

Position 

[455,  55,  485,  85] 

Floating 

off 

Location 

[188,365,512,604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

yonly 

TimeRange 

"auto" 

YMin 

»-5» 

YMax 

Hen 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 

} 

"0" 

Block  { 
BlockType 

Scope 

Name 

"rollrate" 

Position 

[545,  195,  575,  225] 

Floating 

off 

Location 

[188,235,609,474] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-0.05" 

YMax 

"0.01" 

SaveToWorkspace 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

it  i  M 

Samplelnput 

off 

SampleTime 

} 
Block  { 

"0" 

BlockType 

Scope 

Name 

"theta" 

Position 

[475,  100,  505,  130] 

Floating 

off 

Location 

[188,365,512,604] 

Open 

off 
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Grid 

on 

TickLabels 

on 

ZoomMode 

yonly 

TimeRange 

"auto" 

YMin 

".\" 

YMax 

"1.5" 

SaveToWorkspa 

ce 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 
} 

"0" 

Block  { 
BlockType 

Scope 

Name 

"u" 

Position 

[375,  65,  405,  95] 

Orientation 

left 

Floating 

off 

Location 

[188,365,512,604] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"-100" 

YMax 

"0" 

SaveToWorkspa 

ce 

off 

SaveName 

"ScopeData" 

LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation 

"1" 

Samplelnput 

off 

SampleTime 

X 

"0" 

i 
Block  { 

BlockType 

Scope 

Name 

"yawrate" 

Position 

[555,  245,  585,  275] 

Floating 

off 

Location 

[188,235,512,474] 

Open 

off 

Grid 

on 

TickLabels 

on 

ZoomMode 

on 

TimeRange 

"auto" 

YMin 

"0" 

YMax 

"0.05" 

SaveTo  Workspace 

off 

SaveName 

"ScopeData" 
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LimitMaxRows 

on 

MaxRows 

"5000" 

Decimation           "  1 " 

Samplelnput         off 
SampleTime          "0" 

} 
Line  { 

SrcBlock 

"Demux" 

SrcPort 

8 

Points 

[95,  0;  0,  85,  30,  0] 

Branch  { 
Points 
DstBlock 

DstPort 

> 

[0,  -20] 

"To  Workspace8" 

1 

J 
Branch  { 

DstBlock 

"aob" 

DstPort 

} 

1 

} 

Line  { 
SrcBlock 

"Demux" 

SrcPort 

5 

Points 

Branch  { 
Points 
DstBlock 

[50,  0] 

[35,  0;  0,  130] 
"BETADOT" 

DstPort 

} 

Branch  { 
Points 

1 

[0,  -30;  25,  0] 

Branch  { 
Points 

[0,  -25;  40,  0] 

DstBlock 
DstPort 

} 
Branch  { 

DstBlock 

"To  Workspace3 
1 

"beta" 

DstPort 
} 
} 
} 
Line  { 

1 

SrcBlock 

"Demux" 

SrcPort 

4 

Points 
DstBlock 

[40,  0;  0,  -75] 
"theta" 

DstPort 

» 

1 

} 

Line  { 
SrcBlock 

"Demux" 
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SrcPort 

3 

Points 

[0,  -5;  30,  0] 

DstBlock 

"q" 

DstPort 

1 

i 
Line  { 

SrcBlock 

"Demux" 

SrcPort 

2 

Points 

[20,  0] 

DstBlock 

"alfa" 

DstPort 

1 

} 

Line  { 

SrcBlock 

"Demux" 

SrcPort 

1 

Points 

[10,0] 

DstBlock 

"u" 

DstPort 

> 

1 

) 
Line  { 

SrcBlock 

"YAW  FCS" 

SrcPort 

1 

Points 

[5,0] 

Branch  { 

DstBlock 

"Muxl" 

DstPort 

3 

) 
Branch  { 

Points 

[0,  40] 

Branch  { 

DstBlock 

"dr" 

DstPort 
\ 

1 

) 
Branch  { 

Points 

[-50,  0;  0,  20] 

DstBlock 

"To  Workspace  1" 

DstPort 

} 

1 

} 

} 

Line  { 

SrcBlock 

"Mux" 

SrcPort 

1 

DstBlock 

"Ny  " 

DstPort 

1 

} 

Line  { 

SrcBlock 

"Muxl" 

SrcPort 

I 

DstBlock 

"State-space" 

DstPort 

1 
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Line  { 

SrcBlock 

"deltae" 

SrcPort 

1 

Points 

[5,0] 

DstBlock 

"Muxl" 

DstPort 

1 

Line  { 

SrcBlock 

"deltaa" 

SrcPort 

1 

Points 

[0,  50] 

DstBlock 

"Muxl" 

DstPort 

} 

2 

Line  } 
SrcBlock 

"Demux" 

SrcPort 

7 

Points 

[25,  0] 

Branch  { 

Points 

[55,  0] 

Branch  { 

Points 

[-5,  0;  0,  170] 

DstBlock 

"PSI" 

DstPort 

1 

Branch  { 

Points 

[30,  0;  0,  40;  20,  0] 

Branch  { 

DstBlock 

"yawrate" 

DstPort 

1 

Branch  } 

Points 

[30,  0] 

DstBlock 

"To  Workspace5 

DstPort 

1 

} 
} 

Branch  { 
Points 
DstBlock 
DstPort 

} 
Branch  { 

Points 

DstBlock 

DstPort 


[15,  0;0,  155; -290,  0] 
"YAW  FCS" 

3 


[0,  50] 

"RDOT" 

1 


Branch  { 
Points 


[40,  0,  0,  80] 
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DstBlock 
DstPort 


"Mux" 

2 


} 
Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

Branch  { 

DstBlock 
DstPort 

} 

Branch  { 

DstBlock 
DstPort 

} 
} 
Line  { 

SrcBlock 

SrcPort 

Points 

Branch  { 

DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 
} 
Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 
Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

: 


"State-space" 
1 

"Demux" 
1 


"PvDOT" 
1 
[-5,  0] 

"Mux" 

1 


"Scopel" 
1 


"BETADOT" 

1 

[-10,0] 

"Mux" 
3 


[0,  120] 
"Scope" 
1 


"PuIseVnGeneratorl" 

1 

[25,  0] 

"Sum" 

I 


"  Pulse\nGenerator2 " 

1 

[25,  0] 

"Sum" 

2 
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Jne  j 
SrcBlock 

"Ny 

,  » 

SrcPort 

1 

Points 

[-5, 

0] 

Branch  { 

Points 

[-110,  0;  0,-95] 

DstBlock 

"YAW  FCS" 

DstPort 

> 

2 

i 
Branch  { 

Points 

[0,  80] 

Branch  { 

DstBlock 

"Ny" 

DstPort 

1 

} 

Branch  { 

DstBlock 

"To  Workspace2 

DstPort 

1 

Line  { 
SrcBlock 

"Sum" 

SrcPort 

1 

Points 
Branch  { 
Points 

[10,0] 

[0,  40;  -65,  0] 

DstBlock 
DstPort 

} 
Branch  { 

"To  Workspace" 
1 

Points 
Branch  { 

Points 

DstBlock 

[10,  0] 

[5,0] 

"doublet" 

DstPort 

} 
Branch  { 

Points 

DstBlock 

1 

[0,  95] 

"YAW  FCS" 

DstPort 

} 
} 

} 

Line  { 
SrcBlock 

1 

"Demux" 

SrcPort 

6 

Points 
Branch  { 

DstBlock 

[120,  0] 

"rollrate" 

DstPort 

1 
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Branch  { 

Points 

[0,  -20,  25,  0;  0,  -20] 

DstBlock 

"To  Workspace4" 

DstPort 

} 

1 

} 
Line  { 

SrcBlock 

"PSI" 

SrcPort 

1 

Points 

[15, 

0] 

Branch  { 

DstBlock 

"psi" 

DstPort 

1 

Branch  { 

Points 

[0,  -40] 

DstBlock 

"To  Workspace6" 

DstPort 

1 

} 

} 

Line  { 

SrcBlock 

"Clock" 

SrcPort 

1 

DstBlock 

"To 

Workspace?" 

DstPort 

1 

} 
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APPENDIX  C.  MATLAB  CODE  FOR  PLOT  PROGRAM 

plot(time,doublet,'b') 
titleCRudder  Pedal  Deflection') 
xlabel('Time  (sec)') 
ylabel('dpedal  (deg)') 
figure 

plot(time,dr,'-',time,drl,'--',time,dr2,':',time,dr3,'x',time,dr4,'  .',time,dr5,'-.') 

grid 

title(*Rudder  Deflection') 

xlabel('Time  (sec)') 

ylabel('dr  (deg)') 

legend('Baseline','No  Ny  feedback','Shorter  time  constant', 'Linearized  rudder','Reduced  rudder 

gain','Optimum  fix',1) 

figure 

plot(time,Ny,'-',time,Nyl,'-',time,Ny2,':',time,Ny3,'x',tinie,Ny4,'.',time,Ny5,'-.') 

grid 

title('Lateral  Acceleration  at  Commander  Station') 

xlabel(Time  (sec)') 

ylabelCNy  (g)') 

legend('Baseline','No  Ny  feedback','Shorter  time  constant','Linearized  rudder', "Reduced  rudder 

gain', 'Optimum  fix',  1 ) 

figure 

plot(time,beta*  1 80/pi,'-',time,betal  *  1 80/pi,'-- 

',time,beta2*  1 80/pi,':',time,beta3  *  1 80/pi,'x',time,beta4*  1 80/pi,'.',  time,beta5*  1 80/pi,'-.') 

grid 

title('Sideslip') 

xlabel('Time  (sec)') 

ylabelCBeta  (deg)') 

legendCBaseline'^o  Ny  feedback','Shorter  time  constant', 'Linearized  rudder', 'Reduced  rudder 

gain'/Optimum  fix',4) 

figure 

plot(time,r*  1 80/pi,*-',time,r  1  *  1 80/pi,'— 

',time,r2*  1 80/pi,':',time,r3  *  1 80/pi,'x',time,r4*  1 80/pi,'.',time,r5*  1 80/pi,'-.') 

grid 

title('Yaw  Rate') 

xlabel(Time  (sec)') 

ylabel('r  (deg/sec)') 

legendCBaseline'^o  Ny  feedback', 'Shorter  time  constant',*Linearized  rudder'/Reduced  rudder 

gain', 'Optimum  fix',1) 

figure 

plot(time,p*  1 80/pi,'-',time,p  1  *  1 80/pi,'- 

',time,p2*  1 80/pi,':',time,p3  *  1 80/pi,'x',time,p4*  1 80/pi,",time,p5*  1 80/pi,'-.') 

grid 

titleCRoll  Rate') 
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xlabel('Time  (sec)') 

ylabel('p  (deg/sec)') 

legendCBaseline'^o  Ny  feedback'/Shorter  time  constant', 'Linearized  rudder', 'Reduced  rudder 

gain','Optimum  fix',1) 

figure 

plot(time,psi*  1 80/pi,'-',time,psi  1  *  1 80/pi, '-- 

',time,psi2*  1 80/pi,' : ',time,  psi3  *  1 80/pi,  V,time,psi4*  1 80/pi,'.',time,psi5*  1 80/pi,'-.') 

grid 

title('Heading') 

xlabel('Time  (sec)') 

ylabel('psi  (deg)') 

legend(rBaseline',rNo  Ny  feedback'/Shorter  time  constant','Linearized  rudder', 'Reduced  rudder 

gain','Optimum  fix',  1 ) 

figure 

plot(time,aob*  1 80/pi,'-',time,aob  1  *  1 80/pi, '- 

',time,aob2*180/pi,':^time,aob3*180/pi,'x',time,aob4*180/pi,*.•,time,aob5*180/pi,,-.,) 

grid 

titleO  Angle  of  Bank') 

xJabel(Time  (sec)') 

ylabel('phi  (deg)') 

legend(rBaseline',rNo  Ny  feedback'/Shorter  time  constant','Linearized  rudder', deduced  rudder 

gain'/Optimum  fix',4) 
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APPENDIX  D.  STS-62  FLIGHT  DATA  AND  SES  SIMULATION 

RESULTS 

The  first  four  plots  contain  the  STS-62  data.  The  remaining  six  plots  are  the  SES  results. 
For  all  cases,  Main  Gear  Touch  Down  (MGTD)  occurs  at  80.3  seconds  with  derotation  at  91 
seconds.  The  X  scale  reflects  time  from  the  Approach  and  Land  (A/L)  transition,  which  normally 
occurs  around  1 0,000  feet. 
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